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MATHEMATICAL CIRCULATORY SYSTEM 
MODEL 

RELATED APPLICATION DATA 

5 

This application claims the benefit of U.S. Provisional 
Patent Application No. 60/664,723, filed Mar. 24, 2005, 
entitled “Mathematical Circulatory System Model.” This 
application is a continuation-in-part of, U.S. patent applica- 
tion Ser. No. 10/658,638, filed Sep. 9, 2003, now U.S. Pat. No. 10 
7,182,602 entitled “Whole-Body Mathematical Model for 
Simulating Intracranial Pressure Dynamics,” which claims 
the benefit ofU.S. Provisional Patent Application No. 60/409, 
551, filed Sep. 10, 2002, entitled “A Mathematical Model for 
Human Intracranial Pressure and Flow Dynamics,” and U.S. 15 
Provisional Patent Application No. 60/416,508, filed Oct. 7, 
2002, entitled “A Mathematical Model for Human Intracra- 
nial Pressure and Flow Dynamics.” Each of these four prior 
applications are incorporated herein by reference in their 
entirety. 20 

GOVERNMENT LICENSE RIGHTS 

The U.S. Government may have a paid-up license in this 
invention and the right in limited circumstances to require the 25 
patent owner to license others on reasonable terms as pro- 
vided for by the terms of Grant No. DMS-96-26391 awarded 
by the National Science Foundation; and under Grant No. 
NGT5-40110, Grant No. NNG05GH16H, and Cooperative 
Agreement No. NCC5-581 awarded by NASA. 30 

FIELD OF THE INVENTION 

The present invention relates generally to mathematical 
circulatory system models and more particularly relates to a 35 
mathematical circulatory system model including a regula- 
tory mechanism parameter. 

BACKGROUND 

40 

Lumped-parameter models represent an attractive method 
for examining pressure dynamics involving complicated 
human physiology. In a lumped-parameter modeling 
approach, the physiological system is subdivided into a num- 
ber of linked, interacting subunits termed “compartments.” In 45 
general, each compartment contains a single physical con- 
stituent, e.g., blood, cerebrospinalfluid (CSF), or tissue and 
interstitial fluid. However, depending on the model’s com- 
plexity, a given constituent may appear in more than one 
compartment of the model. Dynamics in each compartment is 50 
specified by lumped, time-dependent functions giving com- 
partmental pressures, while incremental changes in flows and 
compartmental volumes are obtained by associating resis- 
tance and compliance parameters with adjacent compart- 
ments. In particular, interaction between adjacent subunits is 55 
assumed to take place at the interfaces of the model’s com- 
partments. 

With few exceptions, previous models of this type have 
adopted restrictions known as the “Kellie-Monro Doctrine” 
to reduce complexity. The Kellie-Monro framework consid- 60 
ers the intracranial system to be completely enclosed within 
the intracranial vault, which is assumed to be rigid. A speci- 
fied inflow of blood to the intracranial arteries provides a 
forcing for the system, and outflow from the jugular bulb is 
assumed to instantaneously equate to this inflow. These 65 
restrictions yield a closed system with constant total volume. 
Strictly intracranial models have produced a number of 
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important results that illuminate the mechanisms of intracra- 
nial pressure adjustments in situations involving both normal 
and pathophysiology. However, the ability of these closed- 
system models to incorporate the influence of important 
extracranial factors on intracranial pressure dynamics is 
clearly limited. For example, the important buffering effects 
of the spinal CSF space on intracranial pressure cannot be 
directly included. From a mathematical point of view, the 
constant volume constraint also produces an over determined 
system of equations that requires special handling to avoid 
singular behavior. 

Lumped-parameter compartmental models of the present 
type have a long history, dating to the earliest such model of 
the intracranial system formulated by Monro in 1783. This 
first model was bi-compartmental, considering incompress- 
ible brain matter and blood as its two constituents. In the work 
of Monro’s student Kellie 40 years later, the vascular com- 
partment was further subdivided into arterial and venous 
blood to produce a three-compartment model. Since the pio- 
neering work of Monroe and Kellie, increasingly more com- 
plex models of the intracranial compartment have been pos- 
ited in order to more realistically describe the relationship 
between intracranial pressures and volumes. There has been a 
steady increase in the number of fluid compartments, the 
introduction of a separate cerebrospinal fluid compartment, 
the inclusion of cardiovascular input, and a relaxation of the 
treatment of system constituents as incompressible fluids and 
matter. As noted in prior art studies, the intracranial system 
involves a number of subsystems, which interact through 
complex mechanisms, and the classical piecewise approach, 
which often considers CSF dynamics separately, is not suited 
to studying multiple parameter changes and the effects of 
interconnected subsystems on each other. By contrast, 
lumped-parameter models of the intracranial system are 
capable of including and linking different subsystems, so that 
such interactions can be examined. 

When considering lumped parameter models, it is impor- 
tant to realize that a compartment does not necessarily corre- 
spond to a precise physical location in the body. For example, 
with a single CSF compartment in a model, CSF in the ven- 
tricles cannot be distinguished from CSF in the subarachnoid 
and spinal spaces. This is one of the main limitations of the 
lumped-parameter approach. Additional spatial resolution 
can be realized only by subdividing the physical system into 
a larger number of compartments based on spatial consider- 
ations. For example, distinct ventricular CSF and extra-ven- 
tricular CSF compartments may be included as opposed to a 
single lumped CSF compartment. In principle, the entire 
body could be finely subdivided in this manner into separate 
compartments to provide the desired degree of spatial reso- 
lution. However, clearly this subdivision process cannot be 
carried to an extreme as the resulting system of linked gov- 
erning equations will rapidly become too laige for practical 
analysis and solution. 

Despite their evolving complexity, two common features 
characterize most earlier lumped parameter models for pres- 
sure dynamics in the intracranial system. The first common 
feature is an assumption that all resistance and compliance 
parameters can be represented by constants. This leads to a 
linear system of governing differential equations. The second 
common feature is adoption of the “Kellie-Monro Doctrine,” 
which assumes that the intracranial system can be confined 
within the cranial vault. By requiring that inflow to the intrac- 
ranial arteries equals outflow from the jugular bulb, this 
assumption produces a closed system that conserves total 
intracranial volume. 
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As mentioned above, when the intracranial space is treated 
as a closed volume-conserving system contained within the 
(nearly) rigid cranial vault, important mechanisms for the 
influence of extracranial physiology on intracranial pressure 
dynamics cannot be included in the resulting models. For 
example, the ability of the spinal portion of CSF space to 
buffer fluctuations of intracranial CSF pressures cannot be 
directly introduced under the Kellie-Monro Doctrine. At least 
two prior art mathematical models of intracranial pressure 
dynamics include aspects of extracranial physiology. The first 
model includes a compliance for CSF storage within the 
lumbar channel. The first model contains three compliances, 
four resistances, and involves differential equations based on 
a hydrodynamic model for the physical system and its elec- 
trical circuit equivalent. The first model allows the dynamic 
relationship between cerebral perfusion pressure, intracranial 
pressure, and cerebral blood flow in various states of auto- 
regulation to be studied. Use of the first model in conjunction 
with clinical data has determined which indices can be 
derived using transcranial Doppler ultrasonography and 
which trends of intracranial pressure and blood pressure are 
useful in clinical tests of autoregulatory reserve. However, 
despite not strictly abiding by the Kellie-Monro Doctrine, the 
first model falls far short of being a whole-body model for 
intracranial pressure dynamics. While the first model can 
include direct buffering effects of CSF within the lumbar 
channel, it does not contain other important extracranial 
physiology. For example, the first model fails to include rep- 
resentations for the extracranial physiology of the cardiovas- 
cular system and the tissues and extracellular fluids of both 
the central and lower body, including cardiovascular auto- 
regulation, colloid osmotic pressure, and a lymphatic system. 

The second model uses a single ground compartment to 
represent the portion of the body below the clavicles. The 
second model contains three resistances, including a resis- 
tance between the intracranial veins and the rest-of-body 
compartments. Thus, outflow from the intracranial system 
depends on pres sure differences and need not instantaneously 
equal the specified inflow. The second model also contains 
three compliances, including a compliance between the CSF 
and rest-of-body compartment that represents the ability of 
the spinal portion of CSF space to expand and buffer CSF 
pressures by compressing the large veins that surround the 
spinal theca. Two of the three model compliances are pressure 
dependent. The associated functional forms for the pressure- 
dependent compliances determined in the second model are 
used in the present work. The second model still falls short of 
being a whole-body model for intracranial pressure dynam- 
ics. The second model merely lumps extracranial physiology 
into a single compartment, rather than finely subdividing it 
into a plurality of compartments. In addition, the second 
model fails to include the external environment as an implicit 
compartment. 

SUMMARY OF THE INVENTION 

In one embodiment, a method of modeling a circulatory 
system is provided. The method includes (a) providing a 
circulatory system model including one or more time-depen- 
dent pressure functions, each of the one or more time-depen- 
dent pressure functions representing a portion of the circula- 
tory system; (b) using a logistic function to represent a 
regulatory mechanism parameter, the regulatory mechanism 
parameter representing a regulatory mechanism having an 
impact on circulatory system function; and (c) solving the 
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logistic function and at least one of the one or more time- 
dependent pressure functions to determine a circulatory sys- 
tem value. 

In another embodiment, a computer readable medium con- 
5 taining computer executable instructions implementing a 
method of modeling a circulatory system is provided. The 
instructions include (a) a first set of instructions for providing 
a circulatory system model including one or more time-de- 
pendent pressure functions, each of the one or more time- 
io dependent pressure functions representing a portion of the 
circulatory system; (b) a second set of instructions for using a 
logistic function to represent a regulatory mechanism param- 
eter, the regulatory mechanism parameter representing a 
regulatory mechanism having an impact on circulatory sys- 
15 tern function; and (c) a third set of instructions for solving the 
logistic function and at least one of the one or more time- 
dependent pressure functions to determine a circulatory sys- 
tem value. 

In still another embodiment, a method of modeling a pres- 
20 sure and volume relationship in a compliant vessel is pro- 
vided. The method includes (a) defining a first parameter as a 
change in pressure within the vessel, the change in pressure 
being due to a contraction of smooth muscles of a wall of the 
vessel; and (b) defining a second parameter as an active com- 
25 pliance for the vessel, the active compliance varying with 
internal pressure, external pressure, and the first parameter. 

In yet another embodiment, a computer readable medium 
containing computer executable instructions implementing a 
method of modeling a pressure and volume relationship in a 
30 compliant vessel is provided. The instructions include (a) a 
first set of instructions for defining a first parameter as a 
change in pressure within the vessel, the change in pressure 
being due to a contraction of smooth muscles of a wall of the 
vessel; and (b) a second set of instructions for defining a 
35 second parameter as an active compliance for the vessel, the 
active compliance varying with internal pressure, external 
pressure, and the first parameter. 

BRIEF DESCRIPTION OF THE DRAWINGS 

40 

For the purpose of illustrating the invention, the drawings 
show aspects of one or more embodiments of the invention. 
However, it should be understood that the present invention is 
not limited to the precise arrangements and instrumentalities 
45 shown in the drawings, wherein: 

FIG. 1 illustrates one example of a simplified schematic of 
one embodiment of a multi -compartmental model for deter- 
mining intracranial pressure dynamics; 

FIG. 2 illustrates one example of a block diagram of the 
50 various compartments of a model according to FIG. 1; 

FIG. 3 illustrates one example of an intracranial arteries 
compartment of the model according to FIG. 1; 

FIG. 4 illustrates one example of a chart of interstitial fluid 
55 pressure vs. relative lymph flow; 

FIG. 5 illustrates one example of a chart of percent reduc- 
tion in Q^ c VS. P A ; 

FIG. 6 illustrates one example of a simplified schematic of 
a computer system including one embodiment of the present 
60 invention; 

FIG. 7 A illustrates one example of a chart of pulsatile 
cardiac forcing (Q^) over time; 

FIG. 7B illustrates one example of a chart of the 120/80 
central artery pressure response over time (solid) and the 
65 100/65 intracranial artery pressure response (dashed); 

FIG. 8A illustrates one example of a chart of the predicted 
CSF pressure vs. volume; 
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FIG. 8B illustrates one example of a chart of conductance 
to CSF output according to one embodiment of the present 
invention; 

FIG. 9 illustrates one example of a chart of the percentage 
of Q^ c drop due to a positional change with all SNS reflexes 
activated (solid) and without any SNS reflexes included 
(dashed); 

FIG. 10 illustrates one example of a chart of percentage 
decrease of P 7 (solid) and of (dashed) as a function of 
percentage blood loss due to hemorrhage; 

FIG. 11 illustrates one example of a chart of the predicted 
percentage of Q AC drop vs. percentage of blood loss due to 
hemorrhage with both SNSz reflexes activated (solid) and no 
SNSz reflexes included (dashed); 

FIG. 12 illustrates one example of simplified schematic of 
another embodiment of a circulatory system model; 

FIG. 13 illustrates one example of a chart of a multiplier of 
Qj vs. systemic indicative pressure, P A -P V ; 

FIG. 14 illustrates one example of a chart of an OVP 
cardiac uptake multiplier vs. P K ; 

FIG. 15 illustrates one example of a chart of an ANS 0 
cardiac uptake multiplier vs. P A ; 

FIG. 16 illustrates one example of a chart of an ANS Z 
fluidity multiplier vs. P^; 

FIG. 17 illustrates one example of a chart of an CNS 0 
cardiac uptake multiplier vs. percent Q 7 ; 

FIG. 18 illustrates one example of a chart of an CNS Z 
fluidity multiplier vs. percent Q 7 ; 

FIG. 19 illustrates one example of a chart of a calibrated 
pressure-volume relationship in arteries during a resting state 
(solid) and maximum smooth muscle contraction state 
(dashed); 

FIG. 20 illustrates one example of a chart of a percent 
arterial pressure (dashed) and cardiac output (solid) with 
respect to a percent blood loss due to hemorrhage; 

FIG. 21 illustrates one example of a chart of a predicted 
ability of arterial pressure to recover following various levels 
of blood loss due to hemorrhage; and 

FIG. 22 illustrates one example of a simplified schematic 
of a computer system including one embodiment of the 
present invention; 

DETAILED DESCRIPTION 

The present disclosure is directed to a mathematical circu- 
latory system model. In one aspect, the system and method of 
the disclosure includes a lumped-parameter model utilizing 
one or more differential equations representing pressure 
dynamics. One exemplary embodiment of the present disclo- 
sure is set forth in greater detail below with reference to FIGS. 
1 to 11. This first embodiment is directed to a whole-body 
mathematical model for simulating intracranial pressure 
dynamics. Another exemplary embodiment of the present 
disclosure is set forth in greater detail below with reference to 
FIGS. 12 to 21. This second embodiment is directed to a 
simplified mathematical circulatory system model simulating 
nervous system regulatory mechanisms. 

Whole-Body Embodiment for Simulating Intracranial Pres- 
sure Dynamics 

In one embodiment a whole-body mathematical model for 
simulating intracranial pressure dynamics is provided. As 
described in U.S. Provisional Patent Ser. No. 60/409,551, 
which is incorporated by reference as if disclosed herein in its 
entirety, the embodiment revokes the Kellie-Monro Doctrine 
by consistently embedding the intracranial system within 
whole-body physiology. In one example of the embodiment, 
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a model includes 17 interacting compartments, of which nine 
lie entirely outside of the intracranial vault. Two distinct 
compartments are defined to distinguish ventricular from 
extraventricular CSF. The vasculature of the intracranial sys- 
5 tern within the cranial vault is also subdivided into five com- 
partments representing fluid in the intracranial arteries, cap- 
illaries, choroid plexus, veins, and venous sinus. The body’s 
extracranial systemic vasculature is divided into six compart- 
ments representing the arteries, capillaries, and veins of the 
to central body and the lower body. Tissue and the associated 
interstitial fluid in the intracranial and lower regions are 
divided into two compartments. A composite compartment 
involving the tissues, organs, and pulmonary circulation of 
the central body and an implicit compartment representing 
1 5 the external environment complete the model . Since the time- 
dependent compartmental pressure functions are obtained 
from physical pressures through a “lumping” procedure that 
involves space-averaging over the physical extent of a com- 
partment, the subdivision of physical constituents into dis- 
20 tinct spatially limited compartments is necessary to provide 
spatial resolution in this modeling approach. 

In addition to allowing direct flows (e.g. arteries to capil- 
laries to veins) between adjacent compartments, the disclo- 
sure includes the extracranial transfer of fluid between capil- 
25 laries and tissue by filtration. An extracranial lymphatic 
system is also included in the model embodiment. Compo- 
nents of the model allow regulation of systemic vascular 
pressures by the sympathetic nervous system, and, at less than 
extreme (high or low) pressures, autoregulation mechanisms 
30 provide constant blood flow to the cerebrovascular capillaries 
and the choroid plexus as well as constant production of CSF 
in the choroid plexus. Fluid intake, renal output of fluid, and 
adjustment of body volume in response to changes in ambient 
environmental pressure are allowed. A realistic representa- 
35 tion for cardiac uptake and cardiac output provides the forc- 
ing for this system. 

The form of the present embodiment is a system of gov- 
erning differential equations for the hilly time-dependent 
compartmental pressure functions. In one aspect, appropriate 
40 forms for the non-constant resistance and compliance param- 
eters in the model, which may be functions of both pressures 
and time, are determined. Calibration of physically realistic 
scale values for parameters and flows is also a step in the 
simulation process of this example. 

45 Turning to the drawings, wherein like reference numerals 
refer to like elements, FIG. 1 illustrates one exemplary 
embodiment of a mathematical model 10 according to the 
present invention. In mathematical model 10, the human body 
is divided into 16 distinct compartments, lettered A, B, C, D, 
50 F, G, I, J, O, P, S, T, V, X, Y, and Z. Compartment A represents 
the intracranial arteries. Compartment B represents the brain. 
Compartment C represents the intracranial capillaries. Com- 
partment D represents the lower capillaries. Compartment F 
represents the ventricular CSF. Compartment G represents 
55 the lower tissue. Compartment I represents the central arter- 
ies. Compartment J represents the central capillaries. Com- 
partment O represents the central veins. Compartment P rep- 
resents the choroid plexus capillaries. Compartment S 
represents the venous sinus jugular veins. Compartment T 
60 represents the extra- ventricular CSF. Compartment V repre- 
sents the intracranial veins. Compartment X represents the 
lower veins. Compartment Y represents the rest of the body. 
Compartment Z represents the lower arteries. Compartment 
M represents the external atmosphere. Compartment H 7 and 
65 H r represent the left and right heart chambers. The physical 
constituents in subunits of the present embodiment are blood, 
CSF, and tissue and interstitial fluid. With a single exception, 
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i.e., compartment Y, each compartment is composed of a 
single constituent. Spatial resolution is obtained by first 
dividing the body into intracranial and extracranial compo- 
nents. To help delineate the intracranial system, a thick line 12 
in FIG. 1 represents a cranial vault 14, which differentiates 
intracranial from extra-cranial compartments. Compartments 
I, Z, Y, J, D, G, O, and X lie completely outside of cranial vault 
14, and compartment T includes both intracranial and extrac- 
ranial physiology. 

As illustrated in FIG. 2, the majority of the compartments 
in the current embodiment are vascular. The 1 1 vascular com- 
partments can be subdivided spatially into three groups: 
intracranial compartments A, C, P, V, and S; central compart- 
ments I, J, and O, and lower compartments Z, D, and X. The 
“lower” compartments represent the region below the pelvis. 
The “central” compartments in the model represent the region 
between the lower body and the clavicles and also include 
extracranial body components above the clavicles but outside 
of the cranial wall. The vascular system in each of these three 
regions, i.e., intracranial, central, and lower, is divided into 
separate artery, capillary, and vein compartments. However, 
in the intracranial space, the choroid plexus capillaries are 
placed in a separate compartment P from the rest of the 
intracranial capillary bed. This allows these capillaries to 
autoregulate to maintain the production of CSF when ven- 
tricular CSF pressures rise or fall. The venus-sinus veins 
where CSF absorption occurs through the arachnoid villa and 
the jugular veins are also placed in a compartment separate 
from the remainder of the intracranial veins, i.e., compart- 
ment S. There are four strictly non- vascular model subunits. 
Two of these compartments represent tissue matter and inter- 
stitial fluid, i.e., lower compartment G and brain compartment 
B, and the other two represent CSF compartments, i.e., ven- 
tricular F and extra-ventricular T. Compartment T is both an 
intracranial and a central compartment. Compartment T, 
which contains extra- ventricular CSF and lies partially out- 
side of the cranial vault, includes both the subarachnoid and 
spinal portions of CSF space. Compartment T serves as a 
bridging compartment and explicitly allows for buffering of 
CSF pressures by the spinal theca. Two of the three regions 
contain a separate compartment for the tissue and interstitial 
fluid. The exception is the central region where the tissue 
matter, organs (except for the heart), interstitial fluid and 
pulmonary circulation are lumped into a composite rest-of- 
body compartment Y. The central region contains an explicit 
representation for the heart pump having a left chamber pump 
H z , and a right chamber pump H^, and a realistic pulsatile 
cardiac output provides the major forcing in the current 
model. External environment compartment M is also an 
implicit 1 7 -th subunit in the model. No attempt has been made 
to depict relative volumes in FIG. 1, and hence the relative 
sizes of these compartments in the figure do not reflect rela- 
tive volumes. 

The pressure dynamics of the intracranial system are gov- 
erned by a system of differential equations within mathemati- 
cal model 10. Four basic assumptions lead to these equations: 

(1 ) all fluids are considered incompressible and isothermal; 

(2) pressure driven flows are laminar and related to pres- 
sure differences by 
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ments i and j respectively, is the lumped resistance, Zy is 
the fluidity (inverse of R^.), the pressure difference P Z> =P.-P 7 . 
and R^—R^.; 

(3) in the case of fluid filtration from the capillaries into the 
5 interstitial space, the flow is governed by the Starling-Landis 

Equation, i.e., 

Filtxation=^ c ,((P -P ^-o ct {% c -%^))=K ct {P c -o c pi ct ) (2) 

where P c is the capillary pressure, P t is the interstitial fluid 
pressure, k c is the blood plasma colloid osmotic pressure, K t is 
the interstitial fluid colloid osmotic pressure, K ct is the filtra- 
tion coefficient, and a ct is the capillary membrane reflection 
coefficient. The notation for pressure difference has been 
1 extended to osmotic pressure differences by defining jz c =jz c - 
7Z t ; and 

(4) the deformation of the membrane between adjacent 
compartments is a function of the change in pressure differ- 
ence between these compartments, i.e., 


dVij ^ d(Pi-Pj) d(P u ) ( 3 ) 

dt ~ ij dt ~ iJ dt 

25 where \C denotes the instantaneous volume of the ‘cup’ 
formed in the membrane at the interface of compartments i 
and j , Cy denotes the compliance between these two compart- 
ments, and C i j=C Ji . 

As follows, the system of governing equations is obtained 
30 by invoking the conservation law 

flow rate in-flow rate out=rate of volume change (4) 

in each compartment. Referring now to FIG. 3, compartment 
35 I from FIG. 1 is illustrated. Compartment I represents the 
central arteries. The flow rate into compartment I is repre- 
sented by Q^jfrom left heart chamber pump H z . The flow rate 
out of compartment I is the sum of the flows associated with 
lumped resistances K IA , R fJ , and R /z , as represented by equa- 
4Q tion ( 1 ) . The rate of volume change of compartment I i s equal 
to the sum of the terms associated with the compliances C/t 
and C IY , as represented by equation (3). 

Specifically, the pressure dynamics of compartment I are 
governed by the following equation: 

45 

dP, T dV Y , ( 5 ) 

Qhi ~ t Z,jP,j +Z lz P lz + Z IA P IA ) = C IT -^---^. 

Note that in equation (5), the subscript “IY” has been 
reversed thereby causing the associated term to be negative. 
50 Also, in both FIGS. 1 and 3, a filled arrow indicates a one-way 
flow and a hollow arrow indicates a pressure dependent resis- 
tance. In terms of pressure differences, the relation of equa- 
tion (4) yields the following additional 12 differential equa- 
tions and three scalar equations for the remaining 1 5 body 
55 compartments: 

Central Capillary Compartment (J): 

Z/jP ij- {KjRP jy~ Vjy^jy) +Z JO P j O )=0; ( 6 ) 

60 Rest of Body Compartment (Y): 


r i~ r j 

Rij ~ Z >j P ‘j’ Q w + Kjy(Pjy - cr jyKjy) - ( Qyo + Qk) 

65 

where Qu is the flow from compartment i into compartment j, 

P z and P ) are the spatially-averaged pressures of compart- 


(7) 


Cym 


dP Y M 

dt 


dVy, 

dt 


+ 


d Vyo 
dt ; 
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Central Venous Compartment (O): 


Extra- Ventricular CSF Compartment (T): 


ZjqPjo + ZxoPxo + ZsoPso + ZjoPto + Qgo + Qyo - Qoh = 


ZftPft + ZbtPbt — ( ZjsPts + ZjoPto ) = 


Lower Artery Compartment (Z): 


Brain Compartment (B): 


ZizPlZ - Zzd P ZD = Czc - 


ZcbPcb +Z fb P fb - ( Z BV P BV +Z BT P BT ) = 


Lower Capillary Compartment (D): 
ZzdP zd~(ZdxP dx+K D g(P dg-Vdg^dg))^ 
Lower Tissue Compartment (G): 


Kdg(Pdg -cdg^dg) - Qgo 


dP GZ dP GX dP GM 

u ZG —5 + <-GX — -7— + ecM — ■ T — , 

dt dt dt 


Lower Venous Compartment (X): 


Zdx Pdx ~ Z x 0 P xo - Cgx ~ 


Intracranial Artery Compartment (A): 


d P a B d Pa F 

Zja P/a ~ ( Z ac Pac + ZapPap) = C ab + Caf—j—\ 

dt dt 


Intracranial Capillary Compartment (C): 
ZacP act(Z C bP cb + ZcvP cf) = ^> 

Choroid Plexus Compartment (P): 


The terms dV r /dt and dV ro /dt in the conservation equa- 
tions for compartments Y, I, and O have been left intact in 
equations (5), (7), and (8). These volume changes include 
20 components that reflect the regulation of arterial pressure by 
the sympathetic nervous system. They are considered further 
below. The terms involving dP^/dt in the conservation equa- 
tions for 6 compartments Y, and G (equations (7) and (11)) 
denote a volume change for which there is no compensation 
25 through an equal but opposite volume change in an adjacent 
body compartment. Specifically, these terms reflect volume 
increases (or decreases) into the ambient environment, which 
is considered infinitely large and unaffected by pressure 
changes within the body. Ambient environment compartment 
30 M is similar to the ground in an electrical circuit model. It will 
be maintained at a constant pressure value in exactly the same 
way that the ground voltage in a circuit model is given a 
prescribed value relative to which all other values are gauged. 
Accordingly, the ambient environmental pressure P w is fixed 
3 5 here at a constant reference value of zero nun Hg . If this is not 
appropriate for a situation being studied, or if the ambient 
pressure changes with time, then V M must be considered as a 
forcing term in the system of equations. 

40 The set of governing equations may be summed to derive 
the following constraint regarding compliances between the 
body and the external environment: 


ZapPap - {Zp F Pp F + Z PV P PV ) = C PF - 


(Qw - Qk ) + ( Qhi - Qoh ) = Gm - 


Intracranial Veins Compartment (V): 


Zcv Pcv + Zpv Ppv + Z B v P B V - ZvsPvS = C B y - 


Venous Sinus — Jugular Veins Compartment (S): 


If cardiac output equals cardiac uptake (Quj = Qoh) an d the 
atmospheric pressure is constant, equation (21) simplifies to 


ZvsPvs + ZjsPts - ZsoPso : 


Ventricular CSF Compartment (F): 


ZppPpp - ( Zp B Pp B + ZpjPpj ) : 


If, in addition, Q ^=0/^=0 or fluid intake equals fluid dis- 
charge (Q^=Q^), then equation (22) implies as expected that 
the net volume change of the entire body system must be zero . 
Furthermore, if Q^>Q^, as will initially be the case when a 
60 glass of water is consumed, the net volume of the body will 
increase. By equation (22), one of the internal compartmental 
pressures: P G , or most likely P r , must now increase. However, 
if as will be the case, C w is large, the resultant pressure 
increase in compartment Y due to the volume input from the 
65 ambient environment will be small. 

The fluidities in equation (1) must be derived, the filtration 
and reflection coefficients in equation (2) must be calibrated, 
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and the compliances in equation (3) must be calibrated before 
model 10’ s governing differential equations may be used in 
simulations. 

Each constant fluidity Z tj between arbitrary adjacent com- 
partments i and j may be calculated by substituting the mean 5 
pressures and the mean flow rates into equation (1) and solv- 
ing for Z„, giving 


3 Pij + p iz 


Q CF , Total capillary filtration into lower tissue 


QpK).35, Total CSF formation 


Qgp=0.001, Imbalance of diffusion on venule side of 
cerebrovasculature 


QFg=0.044, Flow of CSF through the Virkow-Robins 
Space 


Therefore, once the mean flows and pressures are esti- ah n f rpmainina rr 

mated from physical data, the associated constant fluidities ru , u , ® ‘ 

may be calculated from equation (23). Some model fluidities 15 Callbrated based on °^ ua 
are pressure dependent. A discussion of appropriate expres- constant volume assu 

sions for pressure dependent is provided below where CSF, mean state- The call 

cerebrovascular, sympathetic nervous system, and cardiac flows: 
autoregulation mechanisms are modeled. However, even for a Central Rndv Mean FI cm 

non-constant fluidity, a mean scale value may still be calcu- 2Q ^ 

lated from equation (23) . In the calibrations that follow, mean Qoh^Qhi from cardiac < 

flows and pressures reflect physiological values for an aver- 
age human in a lying down (supine) position. Qia=Pu&hi from equati 

To estimate mean flows, it will be assumed that compart- 
mental volumes remain constant in the mean state. Compart- ~Qiz = Pi&m fr° m equation (25) 

mental pressures are certainly pulsatile, so the system will 25 

have a mean state, but not a steady state. However, since Qjj^Qht^ia + Qiz) ^ 0] 

volume changes are related to pressure differences between _ _ 

adjacent compartments, if pressures in adjacent compart- < ^ JO Q IJ ~Q JY ^ om c011: 

ments move in tandem in the mean state, volumes will remain ^ from constant 

unchanged. Therefore, for a given compartment, if all but one 30 

of the mean flows are known, the final mean flow may be Lower Body Mean Flow 
determined from maintaining constant volume. Further, once 

mean flows into and out of a given compartment are deter- Qzd = Q/z from constant 

mined, these values provide data for flows in the adjacent 

compartments. While many flows can be estimated from data 35 Qgo = Qdg from constan 

available in published literature, most mean flow calibrations _ _ _ 

must make use of the constant volume assumption. Qdx~Qzd~Qdg fr° m cc 

As a starting point for the calibration of mean flows, con- _ 

sider the percentages of cardiac output that exit the cen- < ^ Dxfrom constan1 

tral arteries into the three peripheral artery systems: A, I, or Z. 40 T . • i A/r ^ 

_ , E \ n J . 5 5 Intracranial Mean Flows 

1 hese percentages, in decimal torn, are given below. 


All of the remaining mean flows in the model may now be 
calibrated based on equations (24) through (36) by invoking 
the constant volume assumption in each compartment during 
the mean state. The calibration procedure yields the mean 
flows: 

Central Body Mean Flows 

Qoh^Qhi from cardiac output=cardiac input (37) 

Qia =P ia Q H i from equation (24) (3 8) 


Qij-Q H i~(Qia + Qiz) from constant volume in I 


Qjo = Qij- Qjy from constant volume in J 


from constant volume in Y 


Lower Body Mean Flows 

Qzd = Q/z from constant volume in Z 


Qgo = Qdg from constant volume in G 


~Qdx = Qzd~ Qdg from constant volume in D 


Qvo = Qdx from constant volume in X 


=0.15,=% of cardiac output into intracranial region 


p iz =0.25,0.35,=% of cardiac output into lower region 


Pij=l-(p i a +Pi Z ) ,=% of cardiac output into rest of body (26) 

Additional helpful percentages and ratios include 

p p y=0.70,=% of CSF formation from the choroid 

plexus (27) 50 

X=25 0, = Q a c JQ ap (2 8) 

p te =0.80,=% of CSF drained into venous sinus (29) 

Literature values give the following mean flows in ml/min: 55 

Qh/= 5000, 6600, 6900, Cardiac output (30) 


Q Ci ^=2, Total capillary filtration (with 2 A from liver 
and intestines) 


Qap = Q iJ A+ 1 ) from constant volume in A and equa- 
tion (28) 

Qac = ^Qap from equation (28) 

Qpf = 1 } p jQf from equation (27) 

Qcb = { 1 ~Ppf) Q f+Qb v from Qcb~ Qb fK 1 -P p /)Qf 
Qcv^Qac~Qcb from constant volume in C 
Qpv^Qap~Qpf from constant volume in P 
Q vs = Q cv + Qpv + Qb v from constant volume in V 
Qb t = Qcb+Qfb~ Qb v from constant volume in B 
Qft = Qpf~Qfb from constant volume in F 
Qts = P(sQf from equation (29) 


_ _ (2 1 
OjY - Qcf[ 3 + 3 


1 El 

3 Pij Piz 


Total capillary filtration into central tissue 


Q T o = Q-Pts)QF from constant volume in T 


Qso Q vs + Qrs from constant volume in S 


65 Equation (49) states that 70% of the CSF formation comes 
from the choroid plexus and equation (50) states that the 
remaining CSF formation comes as a filtrate from the 10 
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capillaries through the brain via the Virchow-Robins system 
minus the amount reabsorbed at the venule side of the capil- 
laries (Q^ f ). Notice that, as should be the case, Q B r^Q BT =Q r . 

In addition to allowing direct flows (e.g. arteries to capil- 
laries to veins) between adjacent compartments, the present 
embodiment includes the transfer of fluid between capillaries 
and tissue by filtration. These flows are governed by the 
Starling-Landis equation and are driven by differences 
between the colloid osmotic pressures of the blood plasma in 
the capillaries and the interstitial fluid as well as by the usual 
compartmental pressure differences between capillaries and 
tissue. Filtration mechanisms are included between the cap- 
illaries and tissue compartments in the central and lower 
portions of the body in the present embodiment. In the intrac- 
ranial region, significant colloid osmotic pressure differences 
do occur between the intracranial capillary and tissues. How- 
ever, the endothelial cells that make up the intracranial cap- 
illary wall are so tightly joined that not even water molecules 
can usually pass between them. Thus, under normal condi- 
tions, colloid osmotic pressure forces in the intracranial 
region are irrelevant. In the case of highly elevated capillary 
pressure it is possible for the intracranial capillary wall to 
expand enough to allow water molecules to pass between the 
endothelial cells, and at this point both colloid osmotic pres- 
sures and regular pressures start governing filtration rates. At 
these elevated pressures, the description of intracranial cap- 
illary filtration as a function of capillary pressure will clearly 
be nonlinear. However, to simplify the equation for intracra- 
nial filtration, this relationship may still be linearly approxi- 
mated with the slope defined by Z CB , Z BV , and Z FB . 

There are two locations in the embodiment where fluid 
flow is governed by the Starling-Landis equation (2). These 


flows in the central and lower body are 


(Pj-Py) - ajy( Jt j-jt y)) and 

(59) 

Qd g ~Kd g((Pd~P g) ~ Odg&d-^g) ) > 

(60) 


denoting the flow from the capillaries into the tissue subunits 
of the central and lower body regions, respectively. Solving 
these relations for the filtration coefficients in the mean state 
implies 


F _ Qjy 

(( Pj ~ Py) - tr/y(7T/ - n Y )) 

( 61 ) 

^ Qdg 

&DG = — — _ _ 

((Pd ~ Pc ) - v dg(xd ~ n g )) 

( 62 ) 


q^O.8 (65) 

<*zx? =0 - 9 (66) 

The filtration coefficients K JY and in equations (61) and 
(62) are now readily determined. 

Calibration of scale values for all resistance and compli- 
ance parameters from available physical data and other rela- 
tionships must be accomplished before the model’s govem- 
ing differential equations can be used in simulations. A key 
step in calibrating model compliances is determining the 
distensibilities of the extracranial compartments of the 
model. In particular, each extracranial compartment has an 
associated total compliance equal to the product of its asso- 
15 ciated volume and distensibility. It should be noted that cal- 
culating compliances by this technique yields a total central 
artery compliance C 7 = 1.529 ml/mm Hg, which is within 5 
percent of the data value of 1 .445 ml/mm Hg for total arterial 
compliance measured in the ascending aorta by Chemla et al. 
2Q The central venous compliance calculated by these methods 
yields C r =50 ml/nnn Hg while known systemic venous pres- 
sure volume curves suggests a total venous compliance of 50 
ml/mm Hg. Finally, Noordergraaf states that the systemic 
arterial compliance is between 1 and 2 ml/nnn Hg and sys- 
25 temic venous is between 50 and 200 ml/mm Hg. 

In this section, the various compliances in equation (3) that 
relate volume adjustments to pressure differences will be 
calculated. In the intracranial region, compartmental volume 
increases are restricted by the rigid cranial wall. Conse- 
30 quently, compartmental compliances must be pressure 
dependent and diminish to zero as pressure differences 
become large. On the other hand, in extracranial regions, to 
lowest order compliances may consistently be considered 
constant, approximating a linear relationship between pres- 
3 5 sure differences and volume adj ustments . The present intrac - 
ranial pressure-dependent compliances are extensions of 
those derived in the prior art, while the constant extra-cranial 
compliances will be derived from estimations of the volume 
and distensibility of each compartment. 

40 In a simplified four-compartment whole-body model for 
CSF dynamics described in the prior art, there are only two 
pressure difference dependent compliances. They allow vol- 
ume adjustments between the CSF and arterial blood and 
between the CSF and venous blood. These compliances have 
45 the general form described by the relation 

ctj(P u ) = q i e- r u\ p u\ a \ (67) 


The mean flows on the right hand side of these equations 
have been calculated above. Thus, to determine the filtration 
coefficients K JY and K^, it is only necessary to calibrate 
scale values for the mean colloid osmotic pressures and the 
reflection coefficients. 

Mean interstitial fluid colloid osmotic pressure is given by 

ir^%=8 mm Hg (63) 

and blood plasma colloid osmotic pressure by 

n-j=3i D =2 8 mmHg. (64) 

Notice that these values are invariant with respect to central 
or lower regions. This is not the case, however, with the 
reflection coefficients a JY and o DG . The reflection coefficient 
1 1 in the legs is estimated to be approximately 0.9 while the 
coefficient of the upper body is less than this value. This is 
reflected by the assignments 


where P^P^-P p and the subscripts i andj take the values C 7/ 4 
(forarterial/CSF compliance) and C /v 4 (for CSF/venous com- 
pliance). For both pairs of index values, Coefficients and 
parameters in equation (67) are given by 

C^°=6.5333,r /V =0.633431 a^=0.604229 (68) 

C a /=1. 82745, r^yO. 817102 a ay =0. 869393 (69) 

This four-compartment model contains an additional con- 
stant compliance between the CSF and the rest of the body 
(g), which is approximated as 

9^=0. 13333. (70) 

This feature represents the interface of extra-cranial CSF in 
the spinal theca with the rest of the body. It also acts as a 
background compliance so that the total CSF compliance can 
never reach machine zero in simulations. 


50 


55 


60 


65 
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In the present 1 6 compartment model, the division of the 
body’s cerebrospinalfluid (CSF) space is considerably more 
refined. Hence, the three CSF compliances in the simple four 
compartment model of the prior art must be appropriately 
apportioned among the compliances of the present embodi- 
ment. Three decimal percentages may be introduced to 
describe this allocation of C /v 4 : 


p /v =0.164 percentage of C p 4 allocated to C FV 

(71) 

p, o =0.214 percentage of Cy v 4 allocated to C TO 

(72) 

p tt =0.622 percentage of C p 4 allocated to C TS 

(73) 


Equations (71-73) reflect the fact that total CSF volume is 
approximately 140 ml, 23 of which is found in the ventricles, 
30 in the spinal cord subarachnoid space (theca) and the 
remainder in the cerebral cisterns and subarachnoid space. 
Thus, if the distensibility of the membrane is similar in these 
three different components of CSF, then C fv is 23/140=0.164 
of C /v 4 as indicated by p /v in equation (71). This same tech- 
nique yields the values in equations (72) and (73). The com- 
pliance C PF is excluded in this division of Cy v 4 as the choroid 
plexus capillaries are known to dilate and constrict in order to 
maintain a constant pressure difference between these capil- 
laries and ventricular CSF. This maintains a constant genera- 
tion rate of CSF from the choroid plexus at normal pressure 
values. Consequently, a value for the compliance C PF is irrel- 
evant in the absence of a pressure difference change. 

The ratios that lead to equations (71)-(73) also imply the 
following percentages describing the allocation of C a *: 

p a y=0.786 percentage of C a y allocated to C4) 

p z , =0.214 percentage of C a f allocated to C/r C5) 

The background compliance C^ 4 is now divided between 
venous and arterial interfaces based on systemic venous vol- 
ume being four times that of arterial volume, giving 

C^=0.8Cfc 4 (76) 

Carteriar^Cf* (77) 

The CSF -related compliances may now be calculated 
based on the above percentages as: 

ml/mm Hg 

C^(P F ^ 
ml/mm Hg 

Crs(PTS=0-95p„(CA p Ts)+C ve „ o JC T ^1.27626 

ml/mm Hg 

C BI {P BT )=0.05p„(C fi \P BT ) + C vm ^)V BI =0.137057 

ml/mm Hg 

C T o(PTo)=Pto(CA p To)+C vmo J€ T o=0.200936 
ml/mm Hg 

CAAPAF^Pa/CjiPAFl+C^JT^ 0.0261999 
ml/mm Hg 

CiiiP^PilCjiP^+C^JCj^.WSlUll 
ml/mm Hg 

where again, the pressure difference P z — P y is denoted and 

the scale value C l7 =C z;/ .(P .-P^^C^.^.) . The compliance C BT , 
introduced to incorporate volume adjustments between the 
brain and subarachnoid CSF, is allocated 5% of C TS . This 
states that of the bulk intracranial compliance in compartment 
T, 95% is allocated to the interface with the highly distensible 
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venous sinus veins and only 5% to the interface with the less 
distensible brain tissue. A similar allocation is made for the 
bulk intracranial compliance of the ventricular CSF compart- 
ment F with respect to the intracranial veins and the brain 
5 tissue. 

Finally, brain/vascular compliances Cab(Pab) an d C BV 
(P BV ) are defined to similar CSF/vascular compliances as: 

Cab{Pab)=C A f(Pab)C ab = 0.0209523 ml/mm Hg 

to 

C B v(P bv) = C fv^P sr)^!sf^0. 688845 nil/nun Hg. 

When calculating extra-cranial compliances it is advanta- 
geous to first determine the total compliance of each com- 
partment. This quantity will be denoted by C z with the sub- 
15 script indicating the compartment. This type of compliance 
can be described by the relation 


Total < (78) 

20 „ Increase in Comp artmental Volume 

C omp artmental = ; 

Increase in Compartmental Pressure 

Compliance 


25 and may be experimentally determined by measuring the 
pressure change that is induced by an estimated volume 
change and then taking the inverse of this ratio. Unfortu- 
nately, clinical data does not exist for determination of all of 
the compliances required for the current embodiment. How- 
30 ever, there is information in the literature regarding volumes 
and distensibilities. These are related to compliance by 

Compliance=Distensibility- Volume. (7 9) 

Therefore, total compartmental compliance may be calcu- 
35 lated from total compartmental volume (denoted V z ) and 
compartmental distensibility (denoted D z ) through equation 
(79). However, before calculating compartmental volumes 
and distensibilities it is useful to describe how the resulting 
total compartmental compliances will be allocated to the 
40 local inter-compartmental compliances. 

Estimations of the local compliances are made by requiring 
that the sum of the relevant local compliances must equal the 
total compliance being apportioned. With this restriction, the 
embodiment’s predicted compartmental pressure changes 
45 induced by known volume changes should agree with experi - 
mental data. Based on this principle, the apportionment of 
total compliances among local compliances is described by 


Central, Local Compliances: 

50 C /y =Max[Q-C /r ,0] ml/mm Hg (80) 

C ro =Max[C o -C ro ,0] ml/mm Hg (81) 

^KW = Max[Cy-(C / y +C y O ),0] ml/mm Hg (82) 

55 

Lower, Local Compliances: 

Czg=Cz ml/mm Hg (83) 

C ml/mm Hg (84) 

60 

CeM^Max [ C G - (C ZG + C G x) , 0] ml/mm Hg (85) 


Here, the Max[,] function is used to ensure that no negative 
compliances occur in the embodiment. For the particular 
65 choices of total compliances here, this function is not needed, 
but it is included above as it may be required if total tissue 
compliances are calculated to be much smaller. 



US 7,794,230 B2 


17 


Total blood volume (V 77W ) is estimated at 5000 ml and 
5600 ml, where blood plasma makes up about 60% of this. 
The relative volumes of systemic veins, arteries and capillar- 

ies are: 


y sa =0 . 16 V biood systemic artery volume 

(86) 

V sv =0.64V bloo j systemic vein volume 

(87) 

V sc = 0.04V blood systemic capillary volume 

(88) 

y p ubn^D9 V blood pulmonary system volume 

(89) 

Fwu 0 - 07 v blood heart-blood volume 

(90) 


These values are based on estimates that the systemic arter- 
ies comprise 13% of the blood volume while the systemic 15 
arterioles and capillaries comprise 7%. In the present embodi- 
ment, arteriole volume is lumped with the arteries, and hence 
3% of the arteriole/capillary space is shifted here into the 
artery compartment resulting in the percentages above. It 
should also be noted that these percentages imply that sys- 20 
temic venous volume is four times that of systemic artery 
volume. 

As a next step, arterial and venous volumes must be allo- 
cated between the intracranial, lower, and central compart- 
ments. Estimates for this distribution are given below. It 25 
should be emphasized that these are percentages of the sys- 
temic blood and exclude the blood volumes in the heart and 
pulmonary system. 

pV cran =0.07 percentage of systemic blood, intracra- 30 

nial region (91) 

pV Wcr =0.40 percentage of systemic blood, lower 

region (92) 
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Estimations for the pressure-strain modulus are given by 
Milnor and Nichols and O’Rourke for various branches of the 
systemic arteries. From this data, the distensibility of the 
various vessels can be calculated. Of interest here are the 
calculated distensibilities for the ascending aorta and the 
femoral artery, represented in the model by D 7 and D z . Mil- 
nor’s data suggests that D 7 is 0.0036225 mrnHg -1 and D z is 
0.00101 mmHg -1 while Nichols and O’Rourke data suggests 
D 7 is 0.00320141 mmHg" 1 and D z is 0.00237912 mmHg- 1 . 
Averaging these values results in the assignments 

D^O.00341196 (104) 

D^O.001 69456. (105) 

Since it is known that the central veins are eight times as 
distensible as central arteries, D 0 is calculated accordingly by 

D 0 =SD t . (106) 

There are some suggestions in the prior art that the upper 
venous distensibility may be larger that of the lower leg, but 
these results are inconclusive. Therefore, lower venous dis- 
tensibility is set to 

D^D 0 . (107) 

Consistent with prior art findings, the capillary compart- 
ments C, J, and D are considered non-deformable. 

The distensibility of the lower tissue compartment is 
assigned a weighted average of the interfacing vein and artery 
distensibilities. Since the veins have four times the volume of 
the arteries, these weights are used in 


Dz 4" 4Dy 

D g = -=—^ — 1 - 0.0221754 (ml/mm Hg)/ml 


pV ce?I<ra/ =0.53 percentage of systemic blood, central 35 

region (93) 

The above percentages lead to the following volumes for 
the extra-cranial, vascular compartments: 


Vz = P y lower E« = i° wer artery volume 

(94) 40 

Vj=p V cen w U 5a =central artery volume 

(95) 

Vx=p y tower Ev=l° wer vein volume 

(96) 

V 0 =p V centra/ V sv =central vein volume 

(97) 43 

V D =p V Iower U sc =lower capillary volume 

(98) 

Vj=p V centraI V sc =centia.[ capillary volume. 

(99) 


The volumes of the tissue and rest of body compartments 
are now calculated based on interstitial fluid volume and 
intracellular fluid volume. Total interstitial fluid volume is 
estimated in the prior art to be 12 L, while total intracellular 
fluid volume is estimated to be 28 L. Considering that 55 
approximately 2 L of the intracellular fluid is contained in the 
blood cells, the following volume assignments can be made: 

V inter =12000 ml=interstitial fluid volume (100) 

V ifJtra =26000 ml=intracellular fluid volume (101) 60 

With these volume assignments and the percentages in 
equations (92)-(93), estimates for lower tissue and rest-of- 
body volumes become: 

y G =P y low er( Vinter+ E^Hower tissue volume (102) 65 

(103) 


Finally, the distensibility of the rest-of-body compartment 
(Y) is also a weighted average based on its composition: 


PV central^ inter + V intra )D c + (109) 



= 0.0221838 (ml/mm Hg)/ml 


The distensibilities of the extra-cranial compartments have 
now been determined. Combined with the volumes of the 
previous section, each extra-cranial compartment now has an 
associated total compliance equal to the product of its asso- 
ciated volume and distensibility. 

A lymphatic system is also included in the present embodi- 
ment. This system allows interstitial fluid in the tissue com- 
partments of the central and lower regions of the body to flow 
directly into the central venous compartment. Lymphatic 
flows are thus able to balance filtration in these regions and 
establish an equilibrium state. Lymphatic flow is not intro- 
duced in the intracranial region as the brain appears to lack 
lymphatic vessels. However, drainage of interstitial fluid via 
the Virchow-Robins Spaces is accommodated by pathways 
between the Brain and CSF compartments. 

The lymphatic system is represented in the present 
embodiment through the flows Q ro and Q go . These flows 
provide pathways through which interstitial fluid can move 
directly from the interstitium to the central venous system. 
Unlike the flows between adjacent compartments, which are 
driven by pressure differences, the lymphatic flow is gov- 


y Y =P V central V inteA E„*J+ E™/™ =reSt °f body Volume 
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emed almost exclusively by the interstitial fluid pressure. 
This flow increases 20 fold at interstitial fluid pressures near 
4 mm Hg from the corresponding flow at an interstitial pres- 
sure of about -6 nun Hg. Below this interstitial fluid pressure 
value, lymphatic flow ceases. The relative lymphatic curve is 5 
modeled by the logistic expression 

M (110) 

f W* ) = i + ( M -i) e -K»+6) ’ where « = 20 and r = 1 

This function, depicted in FIG. 4 , defines the relative 1C 
lymph flow in terms of interstitial pressure. The large dot in 
FIG. 4 denotes the mean pressure of interstitial fluid in the 
model. Conversion of this relationship to the actual lymphatic 
flows Q go and Q ro is accomplished by defining 


Plymph(Pc)QGO 
QgokV = - 

( 111 ) 

P lymphiP G ) 


n t*\ Ptymph(Py )QyO 
Qyo(t) = _ 

( 112 ) 

P lymph (Py ) 



In this manner, the mean flow is maintained at the mean 
pressure and can increase or decrease accordingly. Based on 25 
equation (110), Q GO can increase from a mean flow of about 
0.24 ml/min to about 4.7 ml/min while Q ro can increase from 
about 1.76 ml/min to about 35 ml/min. 

Blood supply and CSF production is well regulated in the 
human intracranial system. This auto -regulation is achieved 30 
by the dilation and constriction of the blood vessels induced 
by factors such as carbon dioxide concentration, hydrogen 
ion concentration, oxygen concentration and intracranial 
pressures. This embodiment incorporates such auto-regula- 
tion in the form of pressure sensitive fluidities depicting either 35 
vessel constriction in the form of a reduced fluidity or dilation 
in form of increased fluidity. 

In all, the present embodiment contains 16 compliance 
parameters, of which eight are variable and depend on pres- 
sure differences. There are also 23 resistance parameters. 40 
Four intracranial resistances are pressure-dependent and 
autoregulate flow to both the cerebral capillary bed and the 
choroid plexus as well as the production of CSF in the choroid 
plexus. CSF production is known to be nearly constant over a 
wide range of pressures. However, the production of CSF 45 
decreases when perfusion pressure decreases to low values, 
and a mechanism that ramps down CSF production when 
perfusion pressure falls below about 55 mm Hg is also 
included in the model. 

Blood flow into the cerebral capillaries is auto -regulated by 50 
a pressure dependent fluidity %ac defined by 


_ F(P a ) (113) 

* C Pa -P c 55 

where F(P^) is a logistic expression defined by 

60 

M (114) 

F(P a ) = Q ac with M = 1.00001 and 

AC \ + {M- \)e« p A~r a) 

r = 0.25. 

65 

This implies that Q ac = Z ac (^a~^c) = ^ (P^) and cerebral 
blood flow is thus determined almost entirely by intracranial 
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artery pressure and does not diminish significantly until this 
pressure drops below 50 nun Hg. The relationship between 
Q^ c and is displayed in FIG. 5 . 

Blood flow into the choroid plexus is regulated by a pres- 
sure-difference dependent fluidity Z AP defined by 


Zaf ~ (pf-/> P) ' G(/W) 


(115) 


where the multiplier G( P in equation (1 1 5) is a lunction of 

the perfusion pressure ^ per /^ A -^B that is unity when P perf 
exceeds about 55 mm Hg and then falls linearly to zero below 
this value. This implies that 

Q A p=Z AP iP A -P P )=Q AP -G(P perf ). ( 116 ) 

The multiplier G has been included in equation (115) to 
model the fact that CSF production in the choroid plexus 
scales with blood flow, remains nearly constant over a wide 
range of pressures, and is decreased as a function of the 
magnitude of the perfusion pressure when P ' perf falls below 
about 50 mm Hg to about 60 nun Hg. From equation (116), 
blood flow into the choroid plexus remains constant unless 
perfusion pressure falls to a low level. 

Over the wide range of normal pressures, the production of 
CSF is auto -regulated at the venous side of the choroid plexus 
capillaries by the pressure dependent fluidity Z PV defined by 


Zpy = 


QaP ■ G(Pperf) - ZpF ‘ (Pp ~ Pf) 
Pp-Pv 


(117) 


When Pp e ^55, this expression for Z pv will maintain a 
constant pressure difference between the choroid plexus and 
the ventricular CSF. Substituting Z PV with G=1 into the gov- 
erning equation for the choroid plexus, equation (15) reduces 
to 


Cpp- 


dPp F 

dt 


(118) 


Since the compliance C PF must be non- zero to account for 
the known ability of the Choroid Plexus to transmit pressure 
pulsations to the ventricular CSF, the governing equation for 
the choroid plexus compartment becomes simply 


dPpF =Q (119) 

dt 


This implies a constant pressure difference between the 
choroid plexus and ventricular CSF is maintained by equation 
(117) for Pp C ^55 nun Hg. Therefore, 

QpF=Zpp (P p-P p) =Z PF (P P -P F )=Qpp. ( 120 ) 

Since for pressures in the normal range, CSF production in 
the choroid plexus is proportional to constant CSF pro- 
duction from the choroid plexus is thus achieved. Equation 
(119) also eliminates the need to estimate C PF in this model as 
occurrences of this parameter are always multiplied by either 
dP PP /dt or dP^dt. 

Using the above results, it can now be demonstrated how 
Z P F autoregulates CSF production for P per ^55 mm Hg. Sub- 
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stituting Q PF for ZpF •(Pp-P^) in equation (117) and noting 
that Q^p-Qff = Qfv results in the equalities: 


Q PV (12D 5 

Zpv = /w7 

Qpv 

(Pp - Pp ) + (Pp - Py) 

Qpv 10 

Qpp/Zpp + (Pp - Py) 

The last term in this expression reveals the relationship 
between Z PV and P F . Physiologically, Zp V should decrease 
with increasing CSF pressure (P^) causing an increase in 15 
choroid plexus pressure and maintaining constant CSF pro- 
duction across Z PF' It is quite clear from the last equality in 
equation (121) that this is indeed the case. These autoregula- 
tory mechanisms may be substituted into (13) to (1 6) and (1 8) 
to obtain the governing equations for compartments A, C, P, 

V, and F. 20 

The present embodiment also includes a group of regula- 
tory mechanisms associated with the sympathetic nervous 
system (SNS). Two variable resistances in the central and 
lower regions provide for SNS regulation of arterial pressure 
through constriction of the arterioles. Included in these vari- 25 
able resistances is the dramatic “last ditch stand” triggered by 
the SNS when arterial pressure falls to a level where cerebral 
blood flow is affected. The far less extreme SNS regulation of 
arterial pressure through a constriction of the large vascular 
vessels in the central body is also represented in the model. 30 
Active interfaces are placed between the central rest-of-body 
compartment and the central artery and vein compartments. 
When arterial pressure falls, forcing terms in the governing 
equations for compartments Y, I, and O force the volume cups 
at the active Y-I and Y-O interfaces into the vascular compart- 35 
ments, providing the regulatory constriction. An additional 
SNS mechanism in the model regulates central arterial pres- 
sure by increasing the number of heartbeats per minute if 
arterial pressure falls. 

The sympathetic nervous system (SNS) is associated with 40 
reflex mechanisms that act to maintain levels of arterial pres- 
sure and cardiac output if arterial pressure falls. This section 
deals with the portions of the model that represent the capac- 
ity of the SNS to rapidly increase arterial pressure by con- 
stricting both the arterioles and the large vascular vessels. An 45 
embodiment of the regulation of cardiac output by the SNS 
will be developed below. 

A first SNS pressure mechanism involves constriction of 
the large vessels of the circulation in the central body. If 
arterial pressure falls, the large vessels of the circulation in the 50 
central body, especially the veins, strongly constrict to cause 
a rapid increase in arterial pressure. This mechanism is 
included in the current model by placing active interfaces 
between the central compartment Y and the central vascular 
compartments I and O. Forcing terms in the conservation 55 
equation for compartment Y now force the volume cups at the 
interfaces to push into the I and O compartments when arterial 
pressure diminishes, modeling the SNS large vessel constric- 
tion mechanism. 

The conservation equation in compartment Y is described 60 
by 


dVy dVyo dVyi dV YM (122) 

-dr = — + -dr + -dr =Qjy - Qro es 


where the second equality determines the equation. Here, the 
volume change dV YM /d\ between compartment Y and the 
external environment M, as in equation (3), is simply propor- 
tional to the change in pressure difference dP^dt. However, 
dV yo /dt and dV y /dt involve both changes in the pressure 
differences dP YO /dt and dP^dt, respectively, and a forcing 
term describing the active compliance between Y and O and 
between Y and I, respectively. In particular, 


dV YO 

dt 


= Cyo- 


dPn 

dt 


■ + Fyo (0 with Fyo (0 = -24 P' l and 


(123) 


dV Y , 

~dT 


= C Y i 


dP Y , 

dt 


+ F Y i(t) with F Y ,(t) = -8 P', 


(124) 


The governing differential equation in compartment Y is 
now obtained by replacing dV ro /dt and dV J7 /dt in equation 
(7) by the expressions (123) and (124). Similarly, for the 
governing equations in compartments O and I, dV rc /dt is 
replaced in equation (8) by expression (123) and dV r /dt is 
replaced in equation (5) by expression (124), respectively. 

A second SNS pressure regulation mechanism involves 
constriction of the arterioles in most parts of the body (ex- 
cluding the brain and heart) when central artery pressure 
drops. This causes an increase in arterial pressure through an 
increase in the total peripheral resistance. When the arterial 
pressure drop is severe enough to compromise blood flow to 
the brain, this regulatory response is extreme. 

To model this portion of the SNS regulatory response, two 
multipliers are defined for the artery-capillary fluidities Zjj 
and Z^. These variable resistances are of the form 


Zu = =2iL_ . SNSzi (Pi) • SNSz(Qac)- 
Pi-Pj 

and 


(125) 


Zzd = y QzD - • SNSzi (Pi) • SNSz(Qac)- 
Pz - Pd 


(126) 


The first multiplier, 


SNSzi(Pl) = j. 


(127) 


is a function of central body artery pressure that increases 
resistance if P 7 falls below its mean level. The second multi- 
plier, which is a function of the cerebral blood flow Q 4C , is 
defined by 


SNSz{Qac) = 


M 

1 + (M - 1W@a c-Qac) 


where M = 1.1 and 


(128) 


r = 0.02. 


This multiplier models the last-ditch stand when cerebral 
blood flow is significantly reduced. SNSz remains close to 
unity (not activated) until arterial pressure drops to a level 
where Q^ c is affected. At this point, SNSz drops sharply, 
dramatically increasing the resistances R u and P ZD . 

The major source of forcing in the present whole-body 
embodiment comes from the heart and involves the cardiac 
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output and cardiac uptake Q 077 . Because it is known that 
all extra blood that flows into the heart is automatically 
pumped without delay into the aorta and flows again through 
the circulation, the cardiac output will be set equal to cardiac 
uptake, so that QnfQoiz- It I s further known that the normal 5 
heart, functioning without any special stimulation, can pump 
an amount of venous return up to 2.5 times the normal venous 
return before the heart becomes the limiting factor. Therefore, 
a venous return function R is defined that incorporates all of 
the flow into the central venous compartment 10 

R = Qyo+Qgo+Zjo(P j-P o)+Z T o(P t~P o)+Z so {P sr 

Po)+Z X o(PjrP 0 )- (129) 

The cardiac uptake Q 077 is now defined in terms of this 
venous return and two regulatory multipliers SNSo and OVP 15 
by 

Qo^SNSoiP'jyOVPiPoyR (130) 

where P' 7 denotes the time derivative of the arterial pressure 
function P 7 . Since cardiac output is now based on 20 

the venous return through equation (130). 

The multiplier SNSo in equation (130) is associated with 
cardiac regulation by the sympathetic nervous system (SNS). 

It models an SNS mechanism that increases or decreases the 
heart rate in response to a change in arterial pressure. A 25 
decrease in arterial pressure can markedly increase heart 
activity, with the heart rate increasing to as high as 160-180 
bpm from the normal value of 72 bpm. A linear model for the 
SNSo multiplier as a function of the instantaneous time 
derivative P' 7 is 30 


SNSo(P,) = 1 - 


(131) 


35 


This relationship produces an increase in heart rate to about 
twice the normal level as the result of a pressure drop of about 
20 mm Hg. Note that SNSo (0)=1, so that this regulatory 
mechanism is only active when P 7 is changing. 

The OVP function in equation (130) insures that if venous 
pressure drops, then so does cardiac uptake. It is defined by 


ovr\p 0 ) = 


M 

1 + (M - 1 )e KF o- p o) 


where M = 2.5 and r = 0.5. 


(132) 

45 


Since OVP(P 0 )=1 , this regulatory mechanism is not active 50 
when central venous pressure remains at its mean value. 

Using the above cardiac forcing terms leads to a system of 
equations that describe mean pressures in the sense that the 
oscillations that occur about a mean pressure value are 
removed and the mean pressure may be considered a temporal 55 
average over one cardiac cycle. However, when trying to 
resolve circulatory pressure pulsations caused by cardiac out- 
put, a more instantaneous description of cardiac output is 
necessary. Again in this case, cardiac output is set equal to 
cardiac uptake. 60 

Referring now to FIG. 6, it is contemplated that mathemati- 
cal model 10 of FIG. 1 may operate in a computing environ- 
ment including a computer system 20. Computer system 20 
includes a computer central processing unit (CPU) 22, a 
computer memory 24, and input/output devices 26. Math- 65 
ematical model 10 is typically implemented with or embod- 
ied in a computer program 28 which, when executed by 


computing resources within computer system 20, provide the 
functionality of the present invention. Typically, computer 
program 28 resides in computer memory 24 of an individual 
client computer system. Of course, computer program 28 
may reside in the memory of a local or wide area network 
server or in the memory of an equipment computer processor. 
In one embodiment, computer programs 28 include a body 
compartment module 29, a flow equation module 30, and an 
equation solver module 31. Body compartment module 29 
includes instructions for dividing the body into a plurality of 
compartments and a heart pump, each compartment repre- 
senting a portion of the body. Flow equation module 30 
includes instructions for deriving a plurality of differential 
flow equations, each of which correspond to one of the com- 
partments. Equation solver module 31 includes instructions 
for solving the plurality of differential flow equations. The 
instructions within body compartment module 29, flow equa- 
tion module 30, and equation solver module 31 are executed 
within computer programs 28 to simulate the pressure 
dynamics of the intracranial system. Input/output devices 26 
typically include a storage device 36, such as a hard disk 
drive, a keyboard 38, a pointing device 40, i.e., a mouse, and 
a display device 42, such as a monitor. It is contemplated that 
data may be manually input to mathematical model 10 via 
input/output devices 26 or automatically input via patient 
monitoring equipment or the like. 

To validate the present whole-body model embodiment for 
intracranial pressure dynamics, several types of simulations 
assuming normal physiology were carried out and the results 
compared to clinical observations. In the first type of simula- 
tion, a pulsatile form for cardiac output was introduced as 
forcing so that the model’s predicted pressure responses 
could be compared to clinically measured pressure pulsations 
in healthy humans. The pulsatile cardiac output function with 
parameter values n=l 1, (|)=0.27882, and 0^5.958 was used as 
the cardiac forcing function (Q HI ) in the model’s governing 
equations. This function is depicted in FIG. 7 A and represents 
the cardiac output over two cycles. A mean value for central 
artery pressure of 96 mm Hg was prescribed as an initial 
condition for the first simulation. In the second type of simu- 
lation, a constant flow infusion test was used to validate the 
model’s representations for variable intracranial compliances 
and intracranial fluid dynamics. In these simulations, the 
governing equation for the extra-ventricular CSF compart- 
ment (T) was augmented by adding a constant infusion term 
to model the clinical infusion of mock CSF into the lower 
lumbar space. 

With the lower body region explicitly represented by sepa- 
rate compartments in the present embodiment, the effect of 
orthostatic forces on cerebral blood flow as the body changes 
position can be considered. Two additional types of simula- 
tions were now run to validate the modeling of the regulatory 
mechanisms associated with the systemic nervous system. 
The normal physiology value for the resistance R xo was 
increased twofold and R /z was decreased by one half to simu- 
late a change in body position from lying down to standing up, 
and the behavior of the cerebral blood flow Q^ c was then 
examined. These resistance changes were made instanta- 
neously and the modeled effect on cerebral blood flow was 
determined for the next 30 seconds. In the first simulation of 
this type, the model equations included all of the SNS reflexes 
described above. In the second simulation, the SNS terms 
were removed from the governing equations. 

To examine the potential predictive capabilities of the 
embodiment in pathological conditions, simulations were run 
where all cardiac forcing in the model was suddenly stopped 
to simulate cardiac arrest, i.e., and Q 077 were suddenly set 
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equal to zero. A blood volume of about 5600 ml was assumed 
in these simulations, and the response of pressures in the 
embodiment’s circulatory compartments was determined. 
This behavior, and the predicted final circulatory compart- 
mental pressure values, were then compared to clinical results 
associated with cardiac arrest and the filling pressure of the 
circulation. 

As a second example of pathology, simulations of hemor- 
rhagic shock were carried out. Hemorrhage was modeled by 
the inclusion of an outflow path, denoted Q^, from lower 
venous compartment X into ambient atmosphere compart- 
ment M. This flow was calculated so as to achieve an about 
45% loss in blood volume at the end of the simulation. The 
percent changes in central artery pressure, cardiac output, and 
cerebral blood flow were then calculated with respect to per- 
cent blood loss. 

In all of these simulations, the model’s system of differen- 
tial equations was solved numerically using a typical sym- 
bolic mathematical software package employing maximum 
accuracy settings. 

The usual systolic and diastolic values of pulsatile central 
artery pressure are perhaps the best-known values in human 
physiology. FIG. 7B shows the behavior of central artery 
pressure predicted by the present embodiment in response to 
the pulsatile cardiac forcing developed in the prior art. The 
embodiment’ s predicted response to a mean value of 96 mm 
Hg, prescribed as an initial condition, reproduces a peak 
systolic pressure of about 120 mm Hg and a diastolic pressure 
of about 80 mm Hg. FIG. 7B also depicts the intracranial 
artery pressure response. Here, the predicted systolic and 
diastolic pressures are “ 1 00 over 65 .” These values are typical 
of those in the small arteries. 

Because even the most complex mathematical model must 
be based on assumptions and simplifications of actual physi- 
ology, model validation is an essential step in the develop- 
ment process. In the present case, after calibration of param- 
eters associated with healthy human physiology, the model 
was used in two types of simulations and the results compared 
to physical data. In the first type of simulation, the response of 
compartmental pressures to the realistic pulsatile cardiac out- 
put given in FIG. 7A was determined. As shown in FIG. 7B 
the predicted response of central arterial pressure is the typi- 
cal “1 20 over 80” blood pressure readings expected in healthy 
humans. It is worth noting that a mild incisura (dicrotic notch) 
is discemable in this pressure curve. This notch is typical in 
central artery pressure readings, although is it usually more 
pronounced than in FIG. 7B. In particular, the present 
embodiment cannot capture reflected waves which have a 
tendency to cause a second (usually higher) systolic peak and 
hence a more prominent inflection point. 

The response of the intracranial arteries is also shown in 
FIG. 7B . This response agrees well with prior art estimates for 
pressure pulses in the small arteries. The pressure responses 
of other intracranial compartments were also within expected 
ranges. 

Constant flow infusion tests were also simulated using the 
present embodiment. In these experiments, mock CSF was 
infused at a constant rate into the lower lumbar space. The 
pressure of this space was then measured and associated with 
a calculated total CSF volume change. This resulted in deter- 
mination of a curve known as the global pres sure- volume 
relation. The slope of this curve describes the elastance of the 
entire CSF space, including extracranial portions. The inverse 
of the elastance is the more well-know compliance. 

The typical clinical global CSF pressure-volume relation, 
except at extreme pressures, is an S-shaped curve of logistic 
type. It has a lower pressure plateau near resting pressure, 
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defined as the pressure where CSF production in the system is 
just balanced by CSF absorption by the venous system. This 
region of small slope (large compliance) is due to the ability 
of the system to easily accommodate increases in the volume 
5 of the CSF space at these relatively low pres sures through the 
compression of the venous system (excluding the venus-sinus 
veins). As additional CSF volume is added to the system and 
pressures increase, this capacity for adjustment diminishes as 
there will be less venous blood available to eject to further 
l o compres s the veins . Thus , with increasing pres sures , the pres - 
sure-volume curve steepens indicating a reduction in the 
compliance of the system. The slope of the pres sure-volume 
relationship continues to increase for larger infusion volumes 
until the resulting CSF pressures are high enough that the 
1 5 intracranial arteries can begin to be compressed. At this point, 
some additional compliance enters the system. There is a 
point of inflection in the curve followed by a region of 
decreasing slope that leads to an upper pressure plateau at the 
diastolic pressure of the intracranial arteries. Once additional 
20 volume increases beyond this point increase CSF pressures to 
the systolic pressure of the intracranial arteries, there are no 
additional mechanisms available to buffer additional volume 
increases, and the compliance of the CSF system falls to zero. 

The pressure-volume relation obtained by using the 
25 present embodiment to simulate CSF infusion tests is given in 
FIG. 8A. The pressure in this figure is the pressure P r of the 
extraventricular CSF compartment since this compartment 
contains the lower lumber portion of CSF space where clini- 
cal pressure recordings were made. The volume change in 
30 FIG. 8A represents the net change of both ventricular and 
extra- ventricular CSF volumes combined, as this is the quan- 
tity calculated in the experiments. The logistic-like shape of 
the predicted pres sure-volume relationship is maintained 
until CSF pressures reach high values comparable to the 
35 systolic blood pressure. At this point, the curve turns upward 
and compliance rapidly decreases. 

The same infusion simulations that led to the pressure- 
volume relationship in FIG. 8A may also be used to calculate 
a relationship between CSF pressure increases and CSF 
40 absorption. The slope of this relationship is known as the 
conductance of CSF outflow and is denoted by C out . In model 
simulations, this value is easily calculated by 

45 Z re (P r (r)-P s (r))+ (133) 

CZ,W = Z TO (Pr(f) - P 0 (r)) - ( Qts + g ro ) 

P T (t)-Pr 

50 

where the superscript m denotes the model calculation and t 
denotes the time in minutes from the start of the simulated 
infusion. FIG. 5B illustrates the predicted conductance of 
CSF outflow throughout ten minutes of an infusion simula- 
55 tion. 

When C out has been studied experimentally, a linear rela- 
tionship between CSF pressure increase and CSF absorption 
is observed. In prior art experiments, a mean value for C out of 
0.1 1 (ml/min)/mm Hg is given for a sample of eight healthy 
60 volunteers, and it is stated that a value greater than 0.10 is 
probably normal. The values of the conductance of CSF out- 
flow calculated from the present simulations are shown in 
FIG. 8B. These values change with time, but stay within 0.004 
units of the mean value of 0. 1 1 (ml/min)/mm Hg observed in 
65 the prior art. They are also greater than 0.10 as previously 
suggested. Furthermore, the calculated temporal variation of 
Cout™ is sufficiently small that that the relationship between 
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CSF pressure increase and CSF absorption might easily be 
categorized as linear on the basis of clinical data. 

In the current example calibrations, the pressures, flows, 
and hence resistances are determined from data associated 
with the body being in the horizontal (supine) position. Gravi- 
tational influences on the circulation will change when the 
body changes position. Upon standing up, the blood flow into 
the lower body is aided by gravity, but the venous return from 
the lower body is hindered by it. As the lower body in the 
present embodiment is represented by separate compart- 
ments, the effect of a positional change can be considered. To 
simulate the gravity-induced changes associated with stand- 
ing up, the resistance into the lower arteries (R /z ) may be 
decreased by one half while the resistance from the lower 
veins (R TO ) is doubled. This results in pressure and volume 
increases of the lower arteries and veins. In the current simu- 
lations, these resistance changes were made instantaneously 
and the modeled effect on cerebral blood flow was deter- 
mined for the next 30 seconds. In the first simulation, all 
sympathetic nervous system (SNS) reflexes are intact and in 
the second these are removed. As can be seen in FIG. 9, with 
all SNS reflexes activated, cerebral blood flow is maintained 
at 1 00% of its original value. 

Referring now to FIG. 9, the results of modifying the 
resistances between the central and lower body in the model 
to simulate a change in body position from lying down to 
standing up are illustrated. The solid curve in FIG. 9 indicates 
the response of cerebral blood flow to this change in position 
with the SNS reflexes included in the model equations. The 
simulation predicts that with all SNS reflexes activated, cere- 
bral blood flow is maintained at 100% by the model through 
the change in body position. The dashed curve in FIG. 9 
shows the predicted behavior of cerebral blood flow in the 
absence of regulation by the SNS. 

When the SNS reflexes are removed from the model equa- 
tions, cerebral blood flow drops to about 94% of its original 
value due to the positional change. This predicted decrease 
indicates the important role that the SNS regulatory mecha- 
nisms play in the maintenance of cerebral blood flow in the 
model. The decrease also confirms that loss of the SNS 
reflexes can be a factor in orthostatic intolerance. Syncope, or 
fainting, is a common event with many different causes. Syn- 
cope can be due to sudden vasodilatation (vasodepressor or 
“vasovagal” syncope) as well as postural hypotension, in 
which the normal vasoconstrictive reflex response to a tran- 
siently decreased cardiac output on standing is not suffi- 
ciently active. Postural hypotension can occur as the result of 
drugs, venous disease, sympathectomy, hypovolemia, 
peripheral neuropathy, in addition to degeneration of the sym- 
pathetic nervous system (primary autonomic insufficiency, or 
idiopathic orthostatic hypotension). Common to all causes of 
syncope is a decrease in cerebral blood flow to under 30 ml 
per 100 g brain tissue per minute from the usual 50-55 ml, 
about 55% of normal. This is well above the threshold for loss 
of electrical function, at 30% of normal, but is enough to 
cause a transient loss of consciousness. The calculated per- 
centage decrease in cerebral blood flow in the current simu- 
lation with the SNS terms in the model equations omitted 
does not approach the levels associated with fainting due to a 
rapid change in position as the result of standing up too 
quickly, even with inactive SNS reflexes. However, cere- 
brovascular autoregulation in the current simulation remains 
uncompromised and will act to maintain cerebral blood flow 
despite sudden decreases in arterial pressure. 

To demonstrate the potential predictive capabilities of the 
model in pathology, two situations were simulated. In the first 
of these, a simplified representation of cardiac arrest was 
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created by suddenly terminating all cardiac forcing, i.e., set- 
ting Q^Q OH =0, at a specified time in the course of the 
simulation. In prior research by Gutyon, it is noted that 
“When heart pumping is stopped ... the flow of blood every - 
5 where in the circulation ceases a few seconds later. Without 
blood flow, the pressures everywhere in the circulation 
become equal after a minute or so. This equilibrated pressure 
level is called the mean circulatory filling pressure ... at a 
volume 40 of 5000 ml, the filling pressure is the normal value 
to of 7 mm Hg.” Predicted results from the simulation were 
consistent with this statement. When all cardiac forcing was 
stopped in the simulation, all circulatory pressures tend to 
about 7.5 mmHg after about 45 seconds. The small difference 
in the predicted and cited filling pressures may be partially 
15 due to the fact that total blood volume in the simulation was 
taken to be about 5600 ml rather than the 5000 ml. 

Referring now to FIG. 10, the second pathological situation 
considered as a test of the model’s capabilities involves the 
effects of shock caused by hypovolemia. In these simulations, 
20 a flow term Qj M was introduced into the equation for com- 
partment X to model a hemorrhage from the lower body. Prior 
art research by Guyton gives a clinically-derived graph that, 
as in FIG. 10, depicts the percentage drop in central artery 
pressure (P 7 ) and cardiac output (Q^) with respect to percent 
25 blood loss. Consistent with these clinical results, the modeled 
response shows very stable pressure and cardiac output over 
the first 1 0% blood loss, and the relative arterial pressure stays 
above relative output during the entire course of the hemor- 
rhage. At 30% blood loss, the relative drops in pressure and 
30 cardiac output are also in close agreement with the clinical 
results. During the course of the simulated hemorrhage, the 
central venous compartment contracts to less than 50% of its 
original volume, indicating that the active interface (SNSc) in 
the model is causing the central veins to strongly constrict in 
35 order to maintain blood flow back to the heart. Beyond 30% 
blood loss, a second arterial plateau is noticed in FIG. 10 . This 
important feature, which is also noted in the clinical results, is 
due to the “last-ditch stand” reflex to maintain blood supply to 
the brain provided by the regulatory multiplier SNSz (Q 4C ) in 
40 equation (130). 

[Referring now to FIG. 11 , the predicted percentage drop in 
cerebral blood flow (Q^ c ) with respect to percent blood loss 
during the simulated hemorrhage is illustrated. The solid 
curve represents results when both SNSz mechanisms are 
45 intact, and the dashed curve represents results when both 
SNSz mechanisms are disabled. Both simulations predict that 
about 25% blood loss can be tolerated without a significant 
reduction in cerebral blood flow. Above this degree of blood 
loss, with the SNS mechanisms intact, the blood supply to the 
50 brain remains above about 75% of its original value. How- 
ever, without these SNS regulatory mechanisms included in 
the model’s equations, cerebral blood flow drops quickly to 
about 20% its original value. 

The role of the SNS regulatory mechanisms in maintaining 
55 adequate central artery pressure and cardiac output has been 
noted above. However, the need to include the arteriole con- 
strictive reflex SNSz in a model of intracranial pressure 
dynamics when simulating pathology is also dramatically 
illustrated by the behavior of cerebral blood flow as a function 
60 of blood loss in the hemorrhage simulations. Due to the 
intracranial autoregulatory mechanism Z AC in the model, a 
25% blood loss can be tolerated without significant reduction 
in cerebral blood flow. Beyond this point, arteriole constric- 
tion is essential to further maintain adequate levels of cerebral 
65 blood flow. As indicated in FIG. 11 , predicted cerebral blood 
flow Q^ c remains above about 80% of its original value from 
about 25% blood loss through about 40% blood loss. By 
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contrast, in an analogous simulation where the arteriole con- 
strictive response was disabled, Q^ c falls to approximately 
about 55% of its original value when about 40% blood loss 
has occurred. Between about 40% and about 45% blood loss, 

Q 4C remains above about 75% of its original value with the 5 
arteriole constrictive response intact, but falls rapidly from 
about 5 5% to approximately 20% of its original value with the 
arteriole constrictive response disabled. 

Most attempts to study intracranial pressure using lumped- 
parameter models have adopted the classical “Kellie-Monro 10 
Doctrine,” which considers the intracranial space to be a 
closed system that is confined within the nearly -rigid skull, 
conserves mass, and has equal inflow and outflow. The 
present embodiment revokes this Doctrine and develops a 
mathematical model for the dynamics of intracranial pres- 15 
sures, volumes, and flows that embeds the intracranial system 
in extensive whole-body physiology. This new model consis- 
tently introduces compartments representing the tissues and 
vasculature of the extradural portions of the body, including 
both the thoracic region and the lower extremities. In addition 20 
to vascular connections, a spinal subarachnoid cerebrospinal 
fluid (CSF) compartment bridges intracranial and extracra- 
nial physiology allowing explicit buffering of intracranial 
pressure fluctuations by the spinal theca. The embodiment 
may include cerebrovascular autoregulation, regulation of 25 
systemic vascular pressures by the sympathetic nervous sys- 
tem, regulation of CSF production in the choroid plexus, a 
lymphatic system, colloid osmotic pressure effects, and real- 
istic descriptions of cardiac output. Validation tests show that 
the embodiment’s predictions are in agreement with experi- 30 
mental data for normal physiology. Additional simulations 
indicate that the present whole-body model embodiment 
appears to have significant potential predictive capabilities in 
situations involving pathological conditions. 

35 

Simplified Regulatory Mechanisms Embodiment 

In a another embodiment, a logistic function is used in a 
model of a circulatory system to represent a regulatory 
mechanism parameter having an impact on circulatory sys- 
tem function. Such a logistic function may be used in a model 40 
according to the present disclosure. In one example, a logistic 
function according to the present embodiment may be used in 
conjunction with a model such as the whole-body embodi- 
ment described above. In another example, a logistic function 
according to the present embodiment may be used in another 45 
model including one or more time-dependent pressure func- 
tions representing the circulatory system. Solving the logistic 
function and time-dependent pressure function, a desired cir- 
culatory system value may be determined. A modeled circu- 
latory system value can be used in combination with an actual 50 
measured circulatory system value of a subject for a variety of 
applications related to the subject. Example applications 
include, but are not limited to, etiology of Idiopathic Intrac- 
ranial Hypertension, diagnosis of Idiopathic Intracranial 
Hypertension, treatment of Idiopathic Intracranial Hyperten- 55 
sion, modeling behavior of intracranial pressure in micro- 
gravity environments, and any combinations thereof. Further 
discussion of an exemplary microgravity environment appli- 
cation is set forth in U.S. Provisional Patent Application No. 
60/664,723 under the section entitled “Modeling steady-state 60 
intracranial pressures in supine, head-down tilt, and micro- 
gravity conditions.” U.S. Provisional Patent Application No. 
60/664,723 has been incorporated herein by reference in its 
entirety. In one example, a model of the present embodiment 
may be used as an educational tool for generating a model of 65 
a particular physiological system. In such an example a user 
may create a model of a desired physiological system and 
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modify one or more parameters of the model to follow the 
impact that the change would have on the system. In another 
example, a model of the present embodiment may be used in 
conjunction with actual data collected from a subject to be 
studied. Using the real measured data in the model may allow 
prediction of outcomes from the model. These predictions 
may be used to treat the subject. One example of a prediction 
includes the prediction of cerebral blood flow from arterial 
pressure and blood loss data. A simplified example model 
utilizing a logistic function according to the present embodi- 
ment is described in further detail below with respect to FIGS . 
12 to 21. 

In one aspect, new mathematical representations model the 
autonomic and central nervous system reflexes which main- 
tain arterial pressure, cardiac output, and cerebral blood flow. 
In one example, a model may include representations that 
model cardiac up-take, cerebral and non-cerebral blood flow, 
and the pres sure- volume relationship in a vessel with smooth 
muscle contraction. Many of the factors in these relationships 
involve logistic functions. These functions allow regions of 
maximum and minimum effect to be smoothly connected 
through a logistic transition region. 

In another aspect of the present embodiment, a method and 
system of modeling a pressure and volume relationship in a 
compliant vessel is provided. In one example, a first param- 
eter is defined as a change in pressure within a vessel, the 
change in pressure being due to a contraction of smooth 
muscles of a wall of the vessel. A second parameter is defined 
as an active compliance of the vessel, the active compliance 
varying with internal pressure, external pressure, and the first 
parameter. The use of a smooth muscle contraction parameter 
is described in further detail below with a simplified example 
of a model according to the present embodiment. 

Despite its facial simplicity, the present embodiment has 
the ability to accurately capture the physiological responses 
instigated by both the autonomic and central nervous systems 
to regulate cardiac output and arterial pressure and maintain 
the blood flow to the intracranial system necessary to preserve 
essential brain function. Several of these nervous system 
regulatory mechanisms are continuously active and maintain 
a nearly constant arterial pressure under a moderate range of 
stimuli. Others remain inactive until circumstances develop 
that require a triggering of extreme countermeasures to main- 
tain vital blood flow. 

To formulate a lumped-parameter model of a physiological 
system, the system must first be divided into an appropriate 
number of interacting subunits or “compartments.” The num- 
ber of compartments used may depend on a number of con- 
siderations including, but not limited to, the intended appli- 
cation of the model, equation solving resources, the need for 
special resolution, and any combinations thereof. In one 
aspect, a model compartment may not necessarily correspond 
to a precise physical location in the body. For example, in a 
lumped parameter model with a single arterial compartment, 
arterial blood in the thoracic region cannot be distinguished 
from arterial blood in the intracranial arteries. This is one of 
the main limitations of the lumped-parameter approach. 
Because compartmental variables are assumed to be spatially 
averaged over the full extent of a compartment, additional 
spatial resolution can be realized by subdividing the physical 
system into a larger number of compartments. For example, 
to improve spatial resolution, three arterial compartments 
containing blood in the upper, central, and lower body might 
be included as opposed to a single lumped arterial compart- 
ment. In principle, the systemic circulation could be finely 
subdivided in this manner into a sufficient number of separate 
compartments to provide any desired degree of spatial reso- 
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lution. Clearly this subdivision process cannot be carried to 
an extreme as the resulting system of linked governing equa- 
tions will rapidly become too large for practical analysis and 
solution. However, extensive subdivision of the physical sys- 
tem is seldom needed in this modeling approach. A strength 5 
of lumped-parameter models is their ability to represent a 
system over a wide range of parameter values using a reason- 
able (and often small) number of compartments. Further, 
lumped-parameter models not only allow the full extent of the 
physiological system to be studied, but different subsystems 10 
can be consistently linked so that interactions between sub- 
systems, as well as interactions between physical mecha- 
nisms, can be realistically studied. 

Aspects of the present embodiment are described below 
using a simplified model. The example divides the systemic 
circulation into two main compartments representing the 
arteries and the veins. Two flow pathways, representing cere- 
bral blood flow to the brain and non-cerebral blood flow, 
directly connect these two compartments. The example 2Q 
model also contains a simplified heart pump that provides 
forcing for the system. Cardiac output to the arterial compart- 
ment is assumed equal to cardiac up-take from the venous 
compartment. This implies a conservation of blood volume in 
the pulmonary circulation, and allows the model to focus on 25 
the systemic circulation by consistently omitting the pulmo- 
nary system in the model’s description of the cardiovascular 
system. 

Referring to FIG. 12, a schematic 1300 of this example 
model of the present embodiment is illustrated. In this 30 
example, schematic 1300 includes only two main compart- 
ments: an arteries compartment 1305, represented by an A, 
and a veins compartment 1310 represented by the indication 
of a V Each of these compartments has an associated time- 
dependent pressure function, denoted by P A and P^respec- 35 
tively (e.g., measured in millimeters of mercury, “mm Hg”), 
that is spatially averaged over the entire subunit and tempo- 
rally averaged over one cardiac cycle. No attempt has been 
made to depict relative volumes in FIG. 12, and hence the 
relative sizes of these compartments in the figure do not 40 
necessarily reflect relative volumes. Schematic 1300 also 
includes a heart pump 1315 for which cardiac uptake 
(e.g., measured in milliliters per minute, “ml/min”), the inlet 
flow to the right atrium, equals cardiac output Q^, the outlet 
flow from the left ventricle. This constraint imposes a zero 45 
volume change in the pulmonary circulation and allows the 
pulmonary system to be consistently omitted from the exem- 
plary model. Each of arteries compartment 1305 and veins 
compartment 1310 has an associated volume function, 
denoted by V A and V v , respectively. Volumes are related to 50 
pressures and smooth muscle contraction through active com- 
pliance terms c A and C v . The model has two pathways for 
blood to flow directly between the arteries and veins. Cerebral 
blood flow is denoted by Q x and all other flow is denoted by 
Q 2 . This distinction is made because cerebral blood flow, 55 
which may constitute approximately 15 percent of cardiac 
output, is well regulated and remains nearly constant under a 
wide range of arterial pressures to maintain essential brain 
function. While this is also true of blood flow to the heart 
itself, flow to the heart muscle constitutes only about 4 to 5 60 
percent of cardiac output and is separately regulated. The flow 
terms Q^and Q VM are included in the model to allow simu- 
lations of hemorrhage and represent blood flow into a non- 
vascular region. 

Considering blood to be incompressible and balancing 65 
inflows, outflows, and volume adjustments in each compart- 
ment yields the preliminary governing equations 


V a = Qha~Q\-Q?-Qam (134) 

V j^Qx+Qi-Qvm-Qvh (135) 

The dot on the left hand side of these equations designates a 
derivative with respect to time. 

Much of the dynamic behavior described by this embodi- 
ment involves a response that is bounded by a minimum value 
and a maximum value with a smooth logistic transition 
between the two. To aid in the modeling of these responses, 
two preliminary logistic functions are defined in the example 
by 


min -max (136) 

L inc (x, r, max, min) = maxH 7 


min- max (137) 

Ldec (x, r, max, min) = maxH ; : 

1 - min 

1 + e 

max - 1 


where r>0, -oo<x<oo, and min<L zwc ,L rfec <max . The increas- 
ing logistic function L inc tends to max as x — >00 and to min as 
x— >- 00 . The decreasing logistic function behaves in the oppo- 
site manner. Both functions equal unity at x=0, and have an 
inflection point at the average of max and min. 

Cerebral Blood Flow: Q x 

Normal blood flow through the brain is governed by the 
pressure gradient between the cerebral arteries and the brain, 
termed the perfusion pressure. As the present example does 
not have an explicit brain compartment, the pressure differ- 
ence P^-P v , denoted by P AV will act as a surrogate for perfu- 
sion pressure. In this context, P AV will be termed the “sys- 
temic indicative pressure”. Clinical data shows that the flow, 
Q 1? remains nearly constant for arterial pressures between 
about 60 and about 1 50 mm Hg. Casting this in the context of 
flow with respect to systemic indicative pressure, Q x is 
defined here by 

21 =^(^^ 0 . 15 , 1 . 0001 , 0 )-^ ( 138 ) 

where P^is the mean systemic indicative pressure and Q x is 
the mean cerebral blood flow. The multiplier function L inc is 
depicted in FIG. 13. The dot in the figure is located at the 
mean systemic indicative pressure. Q x thus remains nearly 
constant until indicative pressure drops below about 60 mm 
Hg. As venous pressure remains near zero, this threshold 
corresponds to an arterial pressure that is also near 60 m m Hg. 

Non-cerebral Blood Flow: Q 2 

Non-cerebral blood flow, Q 2 , is similarly governed by a 
pressure difference, P AV . Unlike Q l5 it will not be modeled by 
an expression such as in Equation (138). Instead it will be 
governed by the hydrodynamic version of Ohm’s law 

Q 2 =(Pa-Pv)IR2=Z 2 Pav ( 139 ) 

where Z 2 is the inverse of resistance, termed the fluidity or 
conductance. In many mathematical models of this type, flu- 
idities are taken to be constant in a linearization of the gov- 
erning equations. However, in the present example, Z 2 will be 
allowed to vary with both pressures and time. Consequently, 
the embodiment’s governing differential equations (134) and 
(135) will be nonlinear. Z 2 includes three factors; 

Z 2 =ANSz ■ CNSz Z 2 , (140) 

where Z 2 =Q 2 /P AV and the factors ANSz and CNSz vary in 
such a way that they maintain arterial pressure. The functions 
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used to represent ANSz and CNSz are described below in 
relation to details of regulation mechanisms. 

Cardiac Uptake and Output, OVP 

Cardiac output plays a major role in the regulation of 
arterial blood pressure with below-normal arterial pressures 
causing an increase in cardiac output and above-normal pres- 
sures causing a decrease in cardiac output. In the present 
example model, cardiac output has been set equal to 
cardiac uptake Q Venous return and venous pressure are 
two of the major determinants of cardiac uptake. If either is 
allowed to drop significantly, cardiac output will diminish. 

Cardiac uptake is modeled here in terms of the venous 
return Qi+Q 2 -Qj^ and a regulatory multiplier M based on 
central venous pressure and nervous system regulation. In 
particular, 

QvH=MiQl+Q 2 ~QvM)- ( 141 ) 

The dependence of Q VH on the venous return is thus explicit 
in Equation (141), and when M=1 , consistent with the Frank- 
Starling mechanism of the heart. See, Guyton et ah, Textbook 
of Medical Physiology, 10 th Ed. 2000 (cardiac uptake is equal 
to venous return), which is incorporated herein by reference 
in its entirety. The effect of venous pressure on cardiac uptake 
is included in the current model through a factor, OVP (output 
versus pressure), in the cardiac uptake multiplier M. In par- 
ticular, M is defined to be 

M=OVP-ANSo -CNSo ( 142 ) 

where OVP is represented by the logistic function 

OVP(P v )=L inc (Pjr-T v , 0 . 5 , 2 . 5 , 0 ). ( 143 ) 

The behavior of OVP is shown in FIG. 14. The dot in the 
plot is located at the mean venous pressure. Representations 
used in the model for the ANSo and CNSo factors in M are 
described in more detail below. 

Nervous System Regulation of Cardiac Output, Arterial Pres- 
sure, and Cerebral Blood Flow 

Cardiac output and arterial pressure are well regulated. 
Further, since systemic perfusion pressure is most sensitive to 
the relatively large arterial pressure, maintaining arterial 
pressure can insure adequate blood flow to the brain. Two 
levels of regulation are present in the current example of the 
present embodiment. The first level represents the autono- 
mous nervous system (ANS), which provides a regulatory 
response based on changes in arterial blood pressure, trans- 
mitted via the baroreceptors located in the walls of the carotid 
arteries and arch of the aorta. The second level represents the 
ischemic response of the central nervous systems (CNS), 
which is triggered by a significant reduction in cerebral blood 
flow. 

Autonomic Nervous System (ANS) Regulation: ANS effects 
on cardiac output: ANSo 

Changes in the heart rate, and hence the cardiac output, due 
to changes in arterial pressure are included in the present 
example through the factor ANSo in the multiplier M defined 
in equation (142). ANSo thus helps determine how the car- 
diac output differs from the venous return. This factor, 
depicted in FIG. 15, is defined logistically by 

ANSo =L dec (P A -P A , 0 . 1 , 2 , 0 ) ( 144 ) 

where P^ is the normal mean arterial pressure. The dot in the 
plot of FIG. 15 is located at the mean arterial pressure of about 
96 mm Hg. 
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Autonomic Nervous System (ANS) Regulation: ANS Effects 
on Non-Cerebral Blood Flow: ANSz 

ANSz represents the vasoconstriction or vasodilation of 
the arterioles based on changes in arterial pressure. This effect 
5 is included in the example model by varying the fluidity Z 2 in 
equation (140) through the factor ANSz, defined by the logis- 
tic expression 

ANSz=L in JP A -P A , 0 . 3 , 1 . 1 , 0 . 7 ) . ( 145 ) 

10 This function is depicted in FIG. 16. The dot in the plot is 
located at a mean arterial pressure of about 96 mm Hg. 

ANS Effects on Smooth Muscle Contraction in the Walls of 
the Large Vessels 

15 While the two types of ANS reflexes described above will 

aid in maintaining arterial pressure, they will not by them- 
selves necessarily return arterial pressure to its initial mean 
value. Indeed, from a dynamical systems perspective, there 
are infinitely many steady state solutions to equations (134) 
20 and (135) in which arterial pressure is not at its initial mean 
value. An additional ANS regulatory mechanism may there- 
fore be included in the model of the present embodiment to 
ensure that arterial pressure returns to its initial mean value 
under a wide range of stimuli. Smooth muscle contraction in 
25 the walls of the large vessels provides this further mechanism 
by which the ANS regulates arterial pressure and cardiac 
output. This mechanism is incorporated into the current 
example by developing a new form of the pres sure- volume 
relationship in a compliant vessel with smooth muscle con- 
30 traction. 

The total effect of smooth muscle contraction on the pres- 
sure, volume, and compliance of the large vessels is not well 
understood. However, at a constant pressure, smooth muscle 
contraction decreases vessel volume, while smooth muscle 
35 contraction increases pressure at a constant volume. There- 
fore an equation describing the pressure-volume relationship 
in a compliant vessel with smooth muscle contraction may be 
derived from a few basic assumptions. 

For purposes of this example of the embodiment, it is 
40 assumed that each vessel has a maximum volume y Max and a 
minimum volume of zero, and that the transition between 
these extreme values will be logistic in nature. This behavior 
is exhibited by the logistic function relation 

45 

yMax ( 146 ) 


50 where P is the pressure inside the vessel, P* is the pressure 
outside of the vessel, andr>0. The maximum slope, and hence 
maximum compliance, occurs at P=P*. In the present case, 
the ambient pressure P* (e.g., measured in mm Hg) will be 
considered zero. However, P* will be retained in equation 
55 (146) as more complex models may consider these vessels 
residing inside the body where the ambient pressure is 
slightly below zero. 

It is further assumed that the curve defined by V (P) shifts to 
the right during smooth muscle contraction. This is equivalent 
60 to a decrease in volume at a given pressure, or alternatively an 
increase in pressure at a given volume. This shift is a pressure 
increment that will be denoted by (|). Consistent with the 
findings of Bank et ah, “In Vivo Human Brachial Artery 
Elastic Mechanics, Effects of Smooth Muscle Relaxation”, 
65 Circulation 100 (1999) pp. 41-47 and Boutouyrie et ah, 
“Sympathetic Activation Decreases Medium-Sized Arterial 
Compliance in Humans,” Am. J. Physiol 267 (1994) pp. 
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H1368-H1376 (each of which is incorporated herein by ref- 
erence in their entirety) with regard to brachial artery pressure 
and cross-sectional area, it appears that arterial pressure can 
be doubled due to contraction while venous pressure can be 
elevated by as much as about 12 nun Hg. It further appears 
that arterial pressure can be reduced by one half the mean 
arterial pressure and venous pressure can be reduced by about 
8 mm Hg through this mechanism. Therefore, in one 
example, the arteries (|) will be bounded by the range - 
P^/2<(|)^<P^ and by -8<§ v <12 mm Hg in the veins. The 
maximum pressure generated by smooth muscle contraction 
in the arteries will be denoted by ty A Max and in the veins by 
§ v Max . Likewise the minimum pressure generated will be 
denoted by § A min in the arteries and § v min in the veins. It will 
also be assumed that in the resting state, the pressure gener- 
ated by smooth muscle contraction is zero. 

Including the shift (|) due to smooth muscle contraction, the 
expression for V in equation (146) now becomes 


yMax (147) 

1 + e -riP-(<P+P M )) 


where § increases as the smooth muscle contracts. 

Assuming that r and y Max are constants, V can be differ- 
entiated with respect to time to give 


V = C-(P-($ + P*)) 

( 148 ) 

where 


r yMax e -i\P-{4>+P*)) 

( 149 ) 


(1 4 - e - KP-(^+/*))) 2 


and a dot again denotes a time derivative. Here, C represents 
a type of compliance which varies depending on internal 
pressure, external pressure, and the pressure generated from 
smooth muscle contraction. C is referred to as an active com- 
pliance. 

Compliance Calibration 

In compartmental models which do not include smooth 
muscle contraction, volume adjustments are directly related 
to pressure changes through the equation dV/dt=C dP/dt 
where C is a traditional compliance. With smooth muscle 
contraction included, equation (148) shows that this tradi- 
tional relationship will be modified. However, when both the 
pressure shift § and the ambient pressure P* are constants, so 
that (j)=P*=0, equation (148) reduces to the traditional rela- 
tionship and the active compliance C becomes a traditional 
compliance. Consequently, parameter calibration can be 
accomplished in accordance with clinical results where there 
is no change due to smooth muscle contraction. 

The variables in equations (148) and (149) will be identi- 
fied with a subscript indicating the arterial (A) or venous (V) 
compartment. Scale values for x A , r v , V A Max , andV v Max are 
calculated by imposing the following constraints: 


v a (p a $ a ,p* a )=v a , (150) 

Vy(F v j v ,P* V )=V V , (151) 

C a (P a $ a 7 *a)=C a , (152) 

CjAPv,tyv’P*v) = Cv> (153) 
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where, in this model, $ A =§ v =W A =W v =0. The mean volumes 
and compliances on the right hand side of equations (150) to 
(153) are similar to those of the whole body embodiment 
described above. The arterial and venous versions of equation 
5 (148) are used to define the left hand side of the governing 
equations (134) and (135). 

Instead of explicitly defining § A and (|) F , the dynamics gen- 
erated by these terms are included through the introduction of 
two appended differential equations that are solved simulta- 
10 neously with equations (134) and (135). These equations are 

§ A = k A (P A ~P A ) (^ A Max ~^ A )(^. 4~^A mm )( a A^A + ^A) (154) 

^(P A -P A )(ty (0 f”“") ( a (155) 

15 where k A and k v are positive. The last term in each of these 
differential equations is a linear expression in which a and |3 
are determined so that the maximum value of the product of 
the last three terms is unity and this value occurs at the mean 
value <|). Additionally the product of these three terms must 
remain positive for § min ^.§^§ Max . These constraints are sat- 
isfied, when (|)=0, by the following assignments for a and (3 in 
each vessel type: 


_ <p Max + ( f> min (156) 

a ((pMaXfpminf ’ 

p = -((f> Max (t> min y ] - . (157) 

30 

Equations (154) to (157) imply that smooth muscle con- 
traction is most sensitive when (|)=(|). Furthermore, the con- 
traction is initiated by a drop in arterial pressure and is 
bounded between a maximum value of § Max and a minimum 
35 value § mm . Under these conditions, numerical simulations 
imply that appropriate example values for the rate constants 
in equations (154) and (155) are k^O.3 and kj^0.8. 

With the inclusion of equations (154) and (155) above, the 
systemic circulatory system will not come to rest until mean 
40 arterial pressure is achieved or a limiting value of contraction 
is achieved in both types of vessels. Therefore, so long as 
contraction is not at an extremum in both vessels, arterial 
pressure will return to normal. This mechanism can be an 
important factor in maintaining arterial pressure at its mean 
45 value and finalizes the regulation of arterial blood pressure 
via the ANS mechanisms. 

Central Nervous System (CNS) Regulation: Ischemic 
Response 

50 When blood flow through the intracranial region is signifi- 
cantly diminished, the brain can become ischemic, and pow- 
erful responses are rapidly instigated by the central nervous 
system to restore arterial pressure. These responses include 
increased heart rate and vasoconstriction of the arterioles 
55 leading to the non-vital organs. 

CNS Effects on Cardiac Output: CNSo 

The factor CNSo in the multiplier M defined by equation 
(1 42) is associated with a severe change in heart rate triggered 
60 by significant reductions in the cerebral blood flow Q : . This 
effect is included in the model as the determining component 
of cardiac uptake, and is defined logistically by 

CNSo=L dec {Qi-Qi, 0.01,5,0.9) (158) 

65 where Q x is the mean cerebral blood flow. In FIG. 17, CNSo 
is depicted with respect to the percent change in Q x from the 
original value Q x . The dot in the plot indicates a mean value. 
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CNS effects on Non-Cerebral Blood Flow: CNSz 

CNSz is associated with the severe vasoconstriction in the 
arterioles of the non-vital tissues based on changes in cerebral 
blood flow. This effect is included in this embodiment of a 
model as a factor of Z 2 in equation (140) and is defined by 5 

CNSz=L inc (Q i -’Q i ,§.§\,\. §1,0.1). (159) 

In FIG. 18 , CNSz is depicted with respect to the percent 
change in Q x from the original value Q x . The large dot in the 
plot indicates a mean value. 10 

Detailed Governing Equations 

Assuming P*^=P*p?=0, the embodiment’s governing equa- 
tions may now be obtained from the preliminary equations 
(134)-(135), equations (148)-(149), and equations 
(154)-(155). They take the form 15 


Ca'{Pa~^a)~QvH~Qi~Q2~QaM 

(160) 

1^§v) = Qi+Q2-QvM~QvH 

(161) 

$A=k A (PA-P A )(( ^A MaX -^A){^A~^A min )(P-A^A+^A) 

(162) 

^ V = kv(P A~P a) (^l^ aX ~^v)(^J^^V nln )( a J^V + IV) 

(163) 

where 


Q vh =ANSo ■ CNSo OVP ■ (Qi+Q 2 ~Qvm) 

(164) 

Qi=L inc {P A v-PAvJ 0-15, 1.0001, 0)-^ 

(165) 

0 2 = ANSz -CNSz Z 2 - (P A -P V ) 

(166) 

r A V% ax e~ r A (P A- <p A ) 
A (1 + e~ r A( P A-<f , A'>) 2 

(167) 

r v V$ !ax e- r v (P v-4v) 

Cy = t . 

(1 +e~ r v< p v-‘f>v)y 

(168) 


ANSo is defined in equation (144), CNSo is defined in 
equation (158), OVP is defined in equation (143), ANSz is 
defined in equation (145), CNSz is defined in equation (159), 
Q vh is the cardiac uptake, and for the heart pump in the 
present example the cardiac output Qjf A = Q.vir 

Simulations 45 

To test the validity of this example model’s representations 
of the nervous system regulatory mechanisms as well as to 
explore potential applications of the model to pathological 
situations, cardiac arrest and various levels of hypovolemic 
shock were simulated. All simulations were initiated from a 50 
normal resting mean state with values 


P a (0)=P a =96, 

(169) 

Py{ 0)=P^2, 

(170) 

*A<f**S*, 

(171) 

<M0)=<j>^0, 

(172) 


so that the initial value problem involving equations (1 60) to 60 
(163) is completely defined. The value of the cardiac output 
multiplier M and values of the factors ANSz and CNSz in the 
fluidity Z 2 were initially set equal to unity in the simulations 
and then allowed to vary according to equations (142), (145), 
and (159). Values of additional quantities and parameters 65 
used in the simulations are given in Table 1 . The initial value 
problem for the model’s governing nonlinear differential 
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equations was then solved numerically using the mathemati- 
cal software package Mathematica, which is available from 
Wolfram Research Inc., of Champaign, 111. employing maxi- 
mum accuracy and precision settings. 

Cardiac arrest was simulated in the model by setting car- 
diac uptake and output to zero and deactivating all auto- 
regulatory mechanisms. The pressure to which both P 4 and P F 
rapidly tend in this simulation is known as the mean circula- 
tory filling pressure. 

Simulated hemorrhages were studied by assigning a hem- 
orrhage rate to Q^ M and/or Q^In the first case, a total blood 
loss of 45 percent over 30 minutes was simulated, and the 
relative change in arterial pressure and cardiac output over the 
course of the hemorrhage were determined. In a second group 
of simulations, six hemorrhage situations with lesser degrees 
of total blood loss were studied. These cases involved hem- 
orrhaging 24, 28, 31, 36, 41, and 43 percent of total blood 
volume. These simulations were continued beyond the point 
where the hemorrhage was terminated to determine if a recov- 
ery of arterial pressure was predicted with the model’s regu- 
latory mechanisms both active and disabled. Finally, to test 
the ability of the model to accurately capture the overall effect 
of the ANS reflexes, a hemorrhage of 10 percent of blood 
volume was simulated with the ANS mechanisms fully active 
and with the ANS mechanisms deactivated. 


TABLE 1 



Initial and Calibrated Values Used in the Simulations 

Symbol 

Value 

Units 

Description 

?A 

96 

mmHg 

mean (initial) arterial pressure 


2 

mmHg 

mean (initial) venous pressure 

±A 

0 

mmHg 

mean (initial) value of ty A 

<pv 

0 

mmHg 

mean (initial) value of 

Va 

1150 

ml 

mean (initial) arterial volume 

Yn 

4450 

ml 

mean (initial) venous volume 

^VH 

6900 

ml 

mean cardiac uptake 

Q HA 

6900 

ml 

mean cardiac output 

Qi 

0.15 -Q^ 

ml/min 

mean cerebral blood flow 

q 2 

0.85 ■ Q vh 

ml/min 

mean remaining blood flow 

P *A=?*A 

0 

mmHg 

ambient arterial pressure 

* 

III 

'll 

0 

mmHg 

ambient venous pressure 

?A 

2.93657 

ml/mmHg 

mean arterial compliance 

C v 

90.8934 

ml/mmHg 

mean venous compliance 

y A Max 

1318 

ml 

maximum arterial volume 

Vv Max 

8536 

ml 

maximum venous volume 


0.02005 

none 

arterial compliance parameter eq. 
(34) 

*v 

0.04267 

none 

venous compliance parameter eq. 
(35) 

4>A Max 

96 

mmHg 

maximum value of <p A 

$y MaX 

12 

mmHg 

maximum value of 

<1 'A™” 

-48 

mmHg 

minimum value of <p A 

<bv mm 

-8 

mmHg 

minimum value of (|) F 


0.3 

min -1 

arterial contraction parameter eq. 
(21) 

k v 

0.8 

min -1 

venous contraction parameter eq. 
(22) 


Calibrated scale values for the parameters r A , ly, V v Max and 
y^Mcix ca i cu i a t e( i by imposing the constraints (150) to (153) 
are given in Table 1. A portion of the resulting pressure- 
volume relationship in the arteries is depicted in FIG. 19 . As 
a function of arterial pressure V A , the solid curve in FIG. 19 
represents the pressure-volume relationship under resting 
conditions while the dashed curve shows the shift of this 
relationship under maximum smooth muscle contraction. 
The large dot indicates mean pressure and volume in the 
resting state and the arrow indicates the vertical shift associ- 
ated with maximum smooth muscle contraction. 
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Cardiac arrest and mean filling pressure were simulated. 
When cardiac uptake and output were set to zero and all 
auto -regulatory mechanisms deactivated, both P A and P v 
were found to tend towards the same pressure of 6.94 mm Hg 
after about a minute. If the pressures generated by smooth 
muscle contraction and (|) F ) are initially set equal to their 

maximum values and initial arterial and venous pressures are 
adjusted accordingly, the simulations return a value of about 
23 mm Hg, or approximately 3.5 times normal. 

Loss of blood volume by hemorrhage was also simulated. 
The results predicted by simulating a total blood loss of about 
45 percent through hemorrhage with all regulatory mecha- 
nisms active are shown in FIG. 20. This figure depicts the 
relative change in arterial pressure and cardiac output over the 
course of the hemorrhage. Very little difference was noted in 
these simulations when the hemorrhage terms varied between 
arterial hemorrhage and venous hemorrhage. 

Simulations were also conducted with lesser degrees of 
total blood loss to test the ability of the model to predict a 
recovery of arterial pressure once a hemorrhage stops. With 
the model’s regulatory mechanisms active, the predicted 
behavior of arterial pressure during and after hemorrhages 
with total blood losses of 24, 28, 31, 36, 41, and 43 percent of 
total blood volume is shown in FIG. 21. FIG. 21 illustrates the 
predicted ability of arterial pressure to recover following 
various levels of blood loss due to hemorrhage. The resulting 
minimum arterial pressures as a result of these hemorrhages 
were predicted to be approximately 80, 60, 50, 40, 30, and 1 5 
mm Hg, respectively. Recovery of arterial pressure after the 
hemorrhage ceases is predicted in the first three cases, 
whereas no evidence of recovery was seen in the last three 
cases before the simulation was terminated. 

In the simulations of a hemorrhage of 10 percent of total 
blood volume at a rate consistent with FIG. 20, a drop of about 
1 .6 mm Hg in arterial pressure is predicted when all ANS 
mechanisms are active. With all ANS mechanisms deacti- 
vated, the predicted drop is about 61.6 mm Hg. The difference 
between these two cases is thus predicted to be about 60 mm 

Hg. 

The present example is a simplified model having fewer 
compartments than the whole-body model of the embodiment 
described above, which includes sixteen distinct compart- 
ments, a left and right heart, and embeds the intracranial 
system in extensive whole-body physiology. However, 
despite the relative simplicity of the present example, it is 
robust and able to accurately represent the autonomic and 
central nervous system regulatory mechanisms that maintain 
arterial pressure, cardiac output, and cerebral blood flow in 
the face of both moderate and extreme stimuli. In this regard, 
it differs from earlier lumped-parameter models of the car- 
diovascular system. 

Clinical data shows that cerebral blood flow, Q 1; is highly 
regulated and remains nearly constant for arterial pressures 
between about 60 mm Hg and about 1 50 mm Hg. This behav- 
ior is modeled in the present embodiment by a representation 
in which the mean cerebral blood flow, Q l5 is multiplied by a 
logistic function depending on the deviation of the systemic 
indicative pressure P AV from its mean value. As the present 
example has no brain compartment, P^acts as a surrogate for 
perfusion pressure in this expression. As depicted in FIG. 13, 
the model’s representation for Q x remains nearly constant 
until indicative pressure drops below about 60 mm Hg. As 
venous pressure remains near zero, consistent with clinical 
observations, this threshold corresponds to an arterial pres- 
sure that is also near 60 mm Hg. 

As in many lumped-parameter models, the non-cerebral 
blood flow Q 2 is assumed related to the pressure difference 
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^AV through a resistance parameter R 2 , whose inverse is the 
fluidity Z 2 . However, Z 2 here is not assumed to be a constant 
parameter, as is often the case in other lumped-parameter 
models that linearize some or all of the governing differential 
5 equations for the pressure dynamics. Instead, it is a ftmction 
of both pressures and time that involves not only the mean 
scale fluidity value Z 2 but also logistic expressions represent- 
ing the autonomic and central nervous system responses that 
regulate arterial pressure. 

to The factor ANSz in Z 2 has a logistic representation that 
depends on the difference between the arterial pressure P 4 
and its mean value P^. This factor represents the direct regu- 
lation of arterial blood pressure by the autonomous nervous 
system (ANS) via the baroreceptors located in the walls of the 
15 carotid arteries and arch of the aorta. As depicted in FIG. 16, 
ANSz is less than one (vasoconstriction) when arterial pres- 
sure drops below its mean value. This factor now reduces the 
fluidity Z 2 (increases the resistance R 2 ) of non-cerebral flow 
and will cause arterial pressure to increase toward its normal 
20 mean value. Conversely, when arterial pressure ri ses above its 
mean value, ANSz is greater than one (vasodilation), and the 
fluidity Z 2 will be increased, causing arterial pressure to 
decrease toward its mean value. Since parasympathetic 
stimulation does not dilate the arterioles, the only dilation 
25 effect of this mechanism is an inhibition of the vasoconstric- 
tor center. Therefore, the ANS vasodilation has been approxi- 
mated here as one third of the full vasoconstriction capabili- 
ties. Further, the factor ANSz does not diminish to zero as 
arterial pressure drops. Thus, even though essential blood 
30 flow to preserve the integrity of the heart muscle itself has not 
been explicitly included in the non-cerebral flow Q 2 , ample 
flow still remains available for the heart muscle when ANSz 
is at its minimum value. 

The factor CNSz in Z 2 represents part of the ischemic 
35 response of the central nervous system (CNS), which is trig- 
gered by a significant reduction in cerebral blood flow. This 
factor models the severe vasoconstriction in the arterioles of 
the non-vital tissues that is triggered when cerebral blood 
flow drops to levels that can compromise essential brain func- 
40 tions. CNSz is represented by a logistic expression that 
depends on the deviation of Q x from its mean value Q x . As 
shown in FIG. 18, a significant decrease in will induce a 
rapid decrease in the fluidity multiplier CNSz, sharply 
increasing the peripheral resistance of the non-cerebral cir- 
45 culation and promoting an increase of arterial pressure. In 
Guyton et al., this ischemic response has been termed a “last 
ditch stand” in the face of marginal cerebral blood flow. It 
should be noted from FIG. 18 that, as is the case with ANSz, 
the present representation for CNSz also does not diminish to 
50 zero as arterial pressure falls, but tends to a minimum level of 
about 0.1. Thus, even in the face of severely diminished 
arterial pressure, about 4 to 5 percent of the original cardiac 
output will still be available for the heart muscle itself. 

Venous return and venous pressure are two of the major 
55 determinants of cardiac uptake, and hence cardiac output. If 
either is allowed to drop significantly, cardiac output will 
diminish. In the present embodiment, the representation for 
cardiac uptake explicitly involves the venous return. In par- 
ticular, cardiac uptake is defined as the venous return times a 
60 regulatory multiplier M The effect of central venous pressure 
on cardiac uptake enters through this regulatory multiplier, 
which also includes factors that model the regulation of car- 
diac output by the nervous system. 

In the current embodiment, venous pressure affects cardiac 
65 output via a mechanism represented by the factor OVP (out- 
put versus pressure) in the cardiac uptake multiplier. The 
logistic OVP representation models the change in stroke vol- 
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ume associated with a change in central venous pressure. This 
interpretation follows from the fact that when both of the 
nervous system regulatory factors in Mare set to unity, cardiac 
output is not only a multiple of venous return, but the value of 
this multiple is determined entirely by central venous pres- 5 
sure. 

The factor ANSo in M is associated with the autonomic 
nervous system. It models the mechanism that causes the 
heart rate to increase (sympathetic stimulation) when arterial 
pressure drops below its mean value and decrease (parasym- to 
pathetic stimulation) when arterial pressure rises above this 
mean value. The logistic transition portion of ANSo is steep- 
est at the mean value of arterial pressure. Consequently, this 
regulatory mechanism operates most effectively near mean 
arterial pressures. The factor CNSo in the multiplier M is 15 
associated with the ischemic response that produces a severe 
change in heart rate when triggered by significant reductions 
in the cerebral blood flow Ql. The representations for both 
ANSo and CNSo in the present example involve logistic 
functions. 20 

While the autonomic responses modeled by ANSz and 
ANSo will aid in returning arterial pressure toward its normal 
mean value, they may not, by themselves, ensure that arterial 
pressure actually returns to this normal mean value. This is 
because interactions between mechanisms can reduce the 25 
overall effect of the autonomic responses. For example, an 
increase in cardiac activity in response to the ANSo mecha- 
nism will temporarily increase arterial pressure through 
greater cardiac output, but it will also decrease venous pres- 
sure. This, in turn, will reduce the over-all effect of the ANS 30 
response as decreased venous pressure will now cause a drop 
in cardiac uptake and output. Likewise, an increase in sys- 
temic resistance via the ANSz mechanism only ensures a 
greater systemic perfusion pressure. While this will most 
likely entail an increase in arterial pressure it may also result 35 
in a decrease in venous pressure, and again the over-all effect 
will be reduced. An additional ANS regulatory mechanism 
may therefore be included in the model to ensure that arterial 
pressure returns to its initial mean value under a wide range of 
stimuli. 40 

Smooth muscle contraction in the walls of the large vessels 
provides this further mechanism by which the ANS regulates 
arterial pressure and cardiac output. If arterial pressure drops, 
smooth muscle contraction in the arteries tends to return 
arterial pressure back to normal. In addition, contraction in 45 
the large veins increases cardiac output via the OVP mecha- 
nism. This reaction is incorporated into the current embodi- 
ment by defining a new form of the pressure-volume relation- 
ship in a compliant vessel with smooth muscle contraction. 
This relationship gives the effect of smooth muscle contrac- 50 
tion on the pressure, volume, and compliance of the large 
vessels. At a constant pressure, smooth muscle contraction 
decreases vessel volume, while smooth muscle contraction 
increases pressure at a constant volume. Basic assumptions 
now allow compartmental volumes to be represented by a 55 
logistic function involving maximum compartmental vol- 
umes, compartmental and ambient pressures, and a pressure 
increment (|) that reflects the shift of the pressure-volume 
curve to the right during smooth muscle contraction. In this 
representation, maximum compliance is achieved when the 60 
interior pressure P is equal to the sum of the two external 
pressures and P*. 

In the present embodiment, differentiating the pressure- 
volume relationship that results from including the effects of 
smooth muscle contraction produces an expression with an 65 
active compliance parameter that depends on maximum com- 
partment volume, compartmental and ambient pressures, and 
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the contraction pressure increment (|). Unlike prior models, 
the active compliance in this expression multiplies a time 
derivative that includes the ambient pressure P* and pressure 
increment as well as the compartmental pressure P. If the 
ambient pressure is a constant and there is also no change in 
smooth muscle contraction, this expression reduces to the 
usual form. In this case, the active compliance becomes a 
traditional compliance, and this fact simplifies parameter 
calibration. It is interesting to note that the resulting cali- 
brated pressure-volume curves in FIG. 19 modeled according 
to the present embodiment correspond favorably to data col- 
lected by Suga et al. in “Instantaneous Pres sure- Volume Rela- 
tionships and Their Ratio in the Excised, Supported Canine 
Left Ventricle,” Circulation Research 35(1974)pp 117-126, 
which is incorporated herein by reference in its entirety, for 
actual canine left ventricle pres sure- volume relationships. 

To test the ability of the present embodiment to accurately 
represent the autonomic and nervous system regulatory 
mechanisms that maintain arterial pressure, cardiac output, 
and cerebral blood flow, two pathological situations were 
simulated using the example model of FIG. 12 and the results 
compared to clinical observations. In the first of these, cardiac 
arrest was simulated in the model by setting cardiac uptake 
and output to zero and deactivating all regulatory mecha- 
nisms. Within a minute, both V A and P v in the simulation 
tended to the mean circulatory filling pressure. From Rothe, 
“Mean Circulatory Filling Pressure: Its Meaning and Mea- 
surement,”/. Appl. Physiol 74 (1993) pp. 499-509, which is 
incorporated herein by reference in its entirety, the clinical 
value of the normal mean circulatory filling pressure is about 
7 mm Hg. Our simulations return a value of about 6.94 mm 
Hg. If the pressures generated by smooth muscle contraction 
((j)^ and (| ) F ) are initially set equal to their maximum values and 
initial arterial and venous pressures are adjusted accordingly, 
simulations return a value of about 23 mm Hg, or about 3.5 
times normal. This is within the range of about 2 to 4 times 
normal cited in Pang et al., “Peripheral Circulatory Control of 
Preload- Afterload Mismatch with Angiotensin in Dogs,” Br. 
J. Pharmac. 89 (1986) pp. 389-394 and Alexander, “The 
Systemic Circulation,” Annual Rev. ofPhys. 25 (1963) pp. 
213-234, each of which is incorporated herein by reference in 
their entirety. 

The effects of regulation were particularly apparent in 
simulations of the second pathological situation, where hypo- 
volemic shock due to various levels of blood loss through 
hemorrhage was modeled. The first simulation in this set was 
based on a clinical situation involving a 45 percent blood loss 
over the course of 30 minutes. In particular, the modeled 
response displayed in FIG. 20 shows very stable pressure and 
cardiac output over the first 10 percent blood loss, and the 
relative arterial pressure stays above relative output during 
the entire course of the hemorrhage. Beyond about 30 percent 
blood loss, a second arterial plateau is noted in FIG. 20. This 
important feature is due to the triggering of the ischemic 
reflexes that act to maintain blood supply to the brain in the 
face of low perfusion pressure. At about 45 percent blood 
loss, simulation results show cardiac pressure and cardiac 
output dropping quickly to zero. This behavior is in good 
agreement with the measured data on heart rate and mean 
blood pressure given in Opdyke, “Circulatory Effect of Par- 
tial Cerebral Ischemia,” A m. J. Physiol. 143 (1945) pp. 247- 
253, which is incorporated herein by reference in its entirety, 
for the deterioration of the circulation in shock, although a 
small rise in the measured data due to the ischemic response 
prior to a sharp decline occurs in place of the plateau in the 
predicted results. However, the simulation’s predictions, 
including the plateau, for the relative change in arterial pres- 
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sure and cardiac output over the course of the simulated 
hemorrhage with all regulatory mechanisms active are in near 
perfect agreement with the behavior given in Guyton et al. for 
this hemorrhage. 

As noted above, the value of the cardiac uptake multiplier 
Min the simulation of the 45 percent blood volume loss hem- 
orrhage was initially set to unity, consistent with the Frank- 
Starling mechanism of the heart in which cardiac uptake 
equals venous return. The cardiac uptake multiplier was 
found to remain in the range 0.95<M<1.05 over the entire 
course of the simulation, so deviations from the Frank-Star- 
ling “law of the heart” did not exceed about 5 percent. How- 
ever, a separate simulation determined that if M is fixed at 
unity to strictly enforce cardiac uptake equaling venous return 
for all time, a hemorrhage resulting in a blood loss that 
exceeds about 1 2 percent of blood volume will cause death, 
even when all other regulation mechanisms are intact. There- 
fore, while the Frank-Starling “law” is relatively accurate, 
simulations suggest that small deviations apparently must be 
allowed for the maintenance of cardiac output and arterial 
pressure. 

A second group of hemorrhage simulations tested the abil- 
ity of the model to predict recovery from various levels of 
blood loss as well as the critical value of blood loss that 
determines if the hypovolemic shock caused by the hemor- 
rhage will be non-progressive or progressive. Predictions of 
this group of simulations were then compared to the mea- 
sured results presented in Guyton et al., “Dynamics of the 
Heart in Shock,” Fed. Proc. 20, Suppl. 9 (1961) pp. 51-60; 
Walcott, “Blood Volume in Experimental Hemorrhagic 
Shock; "Am. J. Physiol. 143 (1945) pp 247-253; Noble et al., 
“Blood Volume in Clinical Shock: The Extent and Cause of 
Blood of Volume Reduction in Traumatic Hemorrhage and 
Bum Shock,” J. Clin. Invest. 25 (1946) pp. 172-183; and 
Williams et al ..Hematology, 3 rd Ed. (1983), each of which are 
incorporated herein by reference in their entirety. As can be 
seen from FIG. 21 , with all regulatory mechanisms active a 
fairly rapid and complete recovery is predicted by the simu- 
lations for arterial pressures that do not fall below about 50 
mm Hg as a result of the hemorrhage. Conversely, for pres- 
sures that fall to about 40 mm Hg or less, there is little 
recovery made. Since the present embodiment is capable of 
functioning indefinitely at low levels of blood supply to the 
non-vital tissues, progressive shock due to vascular and car- 
diac deterioration is beyond the capabilities of the model. 
However, if arterial pressure remains low in a simulation with 
no evidence of recovery apparent, it is likely that a progres- 
sive shock will occur. Unfortunately, a sharp critical value for 
the minimum arterial pressure level cannot be consistently 
predicted by the current simulations. In the case of minimum 
arterial pressures between about 40 and about 50 mm Hg, 
progression or recovery will most likely depend on long term 
mechanisms such as the renal system, sodium chloride pro- 
duction, and fluid intake which are not represented in the 
model. 

Results from the second group of hemorrhage simulations 
are in close agreement with the measured results in Guyton, 
Walcott, Noble et al., and Williams et al. In Guyton, the “thin 
line of balance between progression of shock and progression 
of recovery” is placed at about 47 nun Hg. The measured data 
in Walcott, Noble et al., and Williams et al. characterizes 
progression or recovery in terms of the percentage of total 
blood volume lost rather than the value of arterial pressure. In 
the experiments of Walcott, recovery following hemorrhage 
occurred for a mean loss of about 34.5 percent of total blood 
volume while an about 43 percent blood loss was progressive. 
In Noble et al., and Williams et al., clinical signs of progres- 
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sive shock occur with a blood loss in the range of about 30 to 
about 40 percent of total blood volume. Consistent with this 
measured data, in the present simulation an about 30 percent 
blood loss results in a minimum arterial pressure of about 54 
5 mm Hg, and from FIG. 21 recovery is predicted. A n about 40 
percent blood loss results in a minimum arterial pressure of 
about 32 mm Hg, and from FIG. 21 progression of shock is 
predicted. 

The importance of the regulatory mechanisms in the mod- 
10 eled hemorrhagic situations can be judged from several fac- 
tors. It is noted in Guyton et al. that with the sympathetic 
reflexes inoperative, recovery does not occur if only about 1 5 
to 20 percent of blood volume is removed over an about 30 
minute period. By contrast, the experimental results Guyton 
15 et al. show that a blood volume loss of between about 30 to 
about 40 percent can be sustained when these reflexes are 
intact and operative. Predictions from the second group of 
hemorrhage simulations agree with these clinical findings. In 
particular, in the simulations, FIG. 21 shows that with the 
20 regulatory reflexes intact, the model predicts that recovery 
can occur after about 30 to about 40 percent blood volume 
loss. On the other hand, when the ANS mechanisms are 
deactivated, the simulations show that recovery is not indi- 
cated at about 15 percent blood volume loss. Of the ANS 
25 mechanisms, the regulation of arterial pressure by smooth 
muscle contraction appears to make the most difference with 
regards to survival following hemorrhage. With smooth 
muscle contraction active but the ANSo and ANSz mecha- 
nisms deactivated, cardiac output and arterial pressure fall to 
30 zero at about a 40 percent blood volume loss (as opposed to 
about 45 percent with all three mechanisms activated). Con- 
versely, if ANSz and ANSo are both active but there is no 
smooth muscle contraction, cardiac output and arterial pres- 
sure fall to zero at about an about 1 8 percent blood volume 
35 loss (as opposed to about 15 percent with all three mecha- 
nisms deactivated). Physiologically, this implies that 
increased heart rate and increased peripheral resistance aid in 
immediate short-term arterial blood pressure maintenance 
but do not act as survival mechanisms. Recovery from a 
40 significant but sub-critical hemorrhage is primarily due to 
regulation of arterial pressure by smooth muscle contraction 
as well as the ischemic reflexes of the central nervous system. 

To insure that the mathematical representations in a model 
of the embodiment accurately capture the overall effect of the 
45 ANS reflex mechanisms, additional simulations involving an 
about 10 percent hemorrhage of total blood volume were 
performed. These simulations predict a difference of about 60 
mm Hg between the case when all ANS mechanisms are 
active and the case when all ANS mechanism are deactivated. 
50 Experiments detailed in Hosomi et al. , “Sinovagal Interaction 
in Arterial Pressure Restoration After 10% Hemorrhage ” Am. 
J. Physiol. Regul. Integr. Comp. Physiol. 237 (1979) pp. 
R203-R209, which is incorporated herein by reference in its 
entirety, determined that following an about 10 percent quick 
55 hemorrhage arterial pressure drops about 7.2 mm Hg with all 
regulatory mechanisms intact and about 67.6 mm Hg with 
blocked reflexes. Therefore, the data shows that these regu- 
latory mechanisms act to reduce the arterial pressure drop 
following the hemorrhage by about 60.4 mm Hg. The model 
60 predictions for this case are thus in excellent agreement with 
measured data indicating that the ANS reflex mechanisms 
have been consistently included in the model and are realis- 
tically represented by their logistic mathematical expres- 
sions. 

65 The methodologies related to this embodiment may be 
implemented using a machine, such as a computer. As dis- 
cussed above with respect to the whole-body model embodi- 
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ment, FIG. 6 illustrates an example computing environment 
for a system and/or method of the present disclosure. It will be 
understood by those skilled in the art that instructions for 
performing methodologies according the present embodi- 
ment may also be implemented in a similar fashion using such 5 
a computing environment. An additional embodiment of a 
machine is provided below with respect to FIG. 22. 

In one aspect, a model of the present disclosure may be 
implemented by a medical device for monitoring a circulatory 
system. In one example, a medical device may utilize the 10 
modeled values of a circulatory system provided by a meth- 
odology of the present disclosure to predict circulatory sys- 
tem behavior. In another example, a medical device may 
utilize modeled values in conjunction with actual measured 
circulatory system values for providing appropriate care to a 15 
patient. A modeled circulatory system value may be com- 
pared to a corresponding circulatory system value that has 
been measured and used in monitoring a circulatory system of 
a subject. 

FIG. 22 shows a diagrammatic representation of one 20 
embodiment of a machine in the exemplary form of a com- 
puter system 2200 within which a set of instructions, for 
causing the machine to perform any one of the methodologies 
of the present invention, may be executed. In alternative 
embodiments, the machine may comprise a network router, a 25 
network switch, a network bridge, Personal Digital Assistant 
(PDA), a cellular telephone, a web appliance, a medical 
device, or any machine capable of executing a sequence of 
instructions that specify actions to be taken by that machine. 
Example medical devices may include, but are not limited to, 30 
an electrocardiograph device, a circulatory pressure monitor- 
ing device, a sphygmomanometer, and any combinations 
thereof. 

The computer system 2200 includes a processor 2205, a 35 
main memory 2210 and a static memory 2215, which com- 
municate with each other via a bus 2220. Computer system 
2200 may further include a video display unit 2225 (e.g., a 
liquid crystal display (LCD) or a cathode ray tube (CRT)). 
Computer system 2200 may also include an alpha-numeric 4Q 
input device 2230 (e.g., a keyboard), a cursor control device 
2235 (e.g., a mouse), a media input device 2240 (e.g., a disk 
drive, a universal serial bus (USB) port, etc.), a signal gen- 
eration device 2245 (e.g., a speaker), and/or a network inter- 
face device 2250. 

45 

Media input device 2240 includes a machine-readable 
medium 2255 on which is stored a set of instructions (i.e., 
software) 2260 embodying any one, or all, of the embodi- 
ments of the present invention. Software 2260 is also shown 
to reside, completely or at least partially, within the main 50 
memory 2210 and/or within the processor 2205. Software 
2260 may further be transmitted or received via the network 
interface device 2250. For the purposes of this specification, 
the term “machine-readable medium” shall be taken to 
include any medium that is capable of storing or encoding a 55 
sequence of instructions for execution by the machine and 
that causes the machine to perform any one of the method- 
ologies of the present invention. The term “machine-readable 
medium” shall accordingly be taken to include, but not be 
limited to, solid-state memories (e.g., random access memory 60 
(RAM), flash memory, etc.), optical and magnetic disks, and 
carrier wave signals. 

Computer system 2200 may also include a circulatory sys- 
tem data import device 2270 for allowing measured data from 
a subj ect to be utilized by set of instructions 2260. Circulatory 65 
system data import device 2270 may be connected to a circu- 
latory system data collection device 2275. 
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Exemplary embodiments have been disclosed above and 
illustrated in the accompanying drawings. It will be under- 
stood by those skilled in the art that various changes, omis- 
sions and additions may be made to that which is specifically 
disclosed herein without departing from the spirit and scope 
of the present invention. 

What is claimed is: 

1. A method of modeling a circulatory system using a 
programmed computer, the method comprising: 

providing a circulatory system model in the computer, the 
model including one or more time-dependent pressure 
functions, each of said one or more time-dependent pres- 
sure functions representing a pressure in a portion of the 
circulatory system; 

using a logistic function in the computer to represent a 
regulatory mechanism parameter, said regulatory 
mechanism parameter representing a regulatory mecha- 
nism having an impact on circulatory system function; 

dividing the circulatory system into a plurality of compart- 
ments and a heart pump, said plurality of compartments 
representing a portion of the circulatory system, said 
heart pump interacting with at least one of said plurality 
of compartments; 

assigning, using the computer, one of said one or more 
time-dependent pressure functions to each of said plu- 
rality of compartments; 

using a smooth muscle contraction parameter, said smooth 
muscle contraction parameter representing an auto- 
nomic nervous system effect on smooth muscle contrac- 
tion in a wall of one of said plurality of compartments, 
wherein said smooth muscle contraction parameter is 
represented according to the following equation: 

V=G(P-(<b+P*)) 

wherein V is a time derivative of a volume of said one of 
said plurality of compartments, P is a time derivative 
of the pressure inside said one of said plurality of 
compartments, $ is a time derivative of said smooth 
muscle contraction parameter, P* is a time derivative 
of the pressure outside of said one of said plurality of 
compartments, and C is an active compliance defined 
according to the following equation: 

r yMax e -riP~{(p + /*)) 

c = T 

(1 4 - e -* P -(^+ P*))) 2 

wherein r>0, y Max is a maximum volume of said one of 
said plurality of compartments, P is the pressure 
inside said one of said plurality of compartments, P* 
is the pressure outside said one of said plurality of 
compartments, and § is said smooth muscle contrac- 
tion parameter; and 

solving, using the computer, one or more equations for said 
one or more time-dependent pressure functions and said 
logistic function to determine a circulatory system 
dynamic. 

2. A method according to claim 1, wherein said dividing 
step divides said circulatory system into a set of compart- 
ments consisting of an arterial compartment and a venous 
compartment, and a heart pump. 

3. A method according to claim 1, wherein said dividing 
step divides said circulatory system into a set of compart- 
ments comprising of an arterial compartment and a venous 
compartment, and a heart pump. 
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4. A method according to claim 3, further comprising using wherein max is a maximum value for x, min is a minimum 
a logistic function to represent a cerebral blood flow param- value for x, r>0, -oo<x<oo, and L^ ec <max. 

eter. 9 . A method according to claim 3, wherein said regulatory 

5. A method according to claim 4, wherein said using of a mechanism parameter represents an autonomic nervous sys- 


logistic function to represent said cerebral blood flow param- 
eter comprises representing said cerebral blood flow param- 
eter according to the following logistic function: 

21 =^(^^ 0 . 15 , 1 . 0001 , 0 )-^ 

wherein Q x is said cerebral blood flow parameter, P "av is a 
pressure difference P A -P V , wherein P A is said time-dependent 
pressure function for said arterial compartment and P K is said 
time-dependent pressure function for said venous compart- 
ment, P^ F is a mean systemic indicative pressure, Q x is a mean 
cerebral blood flow, and L inc is an increasing logistic function 
represented by the following equation: 

min - max 

Line (x, r, max, min) = max H ; : 

1 - min 

1 + -e™ 

max - 1 

wherein max is a maximum value for x, min is a minimum 
value for x, r>0, -oo<x<oo, and min<L Iwc . 

6. A method according to claim 3 , further comprising using 
a logistic function to represent an output versus pressure 
parameter, said output versus pressure parameter represent- 
ing an effect of venous pressure on cardiac uptake. 

7. A method according to claim 6, wherein said using of a 
logistic function to represent an output versus pressure 
parameter comprises representing said output versus pressure 
parameter according to the following logistic function: 

OVP(P v )=L inc (Pjr-P v , 0.5, 2.5,0) 

wherein OVP is said output versus pressure parameter, P K is 
said time-dependent pressure function for said venous com- 
partment, P v is the mean systemic venous pressure, and L inc is 
an increasing logistic function represented by the following 
equation: 


min - max 

L; nc (X r, max, mm) = max H : : 

1 - min 

1 H -e™ 

max - 1 

wherein max is a maximum value for x, min is a minimu m 
value for x, r>0, -oo<x<oo 5 and min<L zwc . 

8. A method according to claim 3, wherein said regulatory 
mechanism parameter represents an autonomic nervous sys- 
tem effect on cardiac output according to the following equa- 
tion: 

ANSo=L dec (P A -P A , 0.1, 2,0) 

wherein ANS 0 is said regulatory mechanism parameter rep- 
resenting said autonomic nervous system effect on cardiac 
output, P A is said time-dependent pressure function for said 
arterial compartment, V A is an initial mean arterial pressure, 
and L dec is a decreasing logistic function represented by the 
following equation: 


min - max 

Ldedx, r, max, min) = maxH : 

1 — min 

1 + rr“ 

max - 1 


5 tern effect on non-cerebral blood flow according to the fol- 
lowing equation: 

ANSz=L inc {P A -P A ,0.3,1 . 1 ,0.7) 

wherein ANSz is said regulatory mechanism parameter rep- 
10 resenting said autonomic nervous system effect on non-cere- 
bral blood flow, P A is said time-dependent pressure function 
for said arterial compartment, is an initial mean arterial 
pressure, and L inc is an increasing logistic function repre- 
sented by the following equation: 

15 

min -max 

LinAx, r, max, min) = maxH ; ; 

mcv ' 1 - min 


wherein max is a maximum value for x, min is a minimum 
value for x, r>0, -oo<x<oo, and min<L zwc . 

1 0. A method according to claim 3, wherein said regulatory 
25 mechanism parameter represents a central nervous system 

effect on cardiac output according to the following equation: 

CNSo=L dec {Q l -'Q 1 , 0 . 01 , 5 , 0 . 9 ) 

wherein Q : is cerebral blood flow, Q x is mean cerebral blood 
30 flow, and L dec is a decreasing logistic function represented by 
the following equation: 

min - max 

Ldec {x, r, max, min) = maxH : 

k 1 — min 

35 1 + -er" 

max - 1 

wherein max is a maximum value for x, min is a minimum 

value for x, r>0, -oo<x<, and L , <max. 

40 . . 

11 . A method according to claim 3, wherein said regulatory 

mechanism parameter represents a central nervous system 
effects on non-cerebral blood flow according to the following 
equation: 

45 CiV5z=i I - MC (e 1 -e 1 ,0.01,1.01,0.1) 

wherein Q x is cerebral blood flow, Q : is mean cerebral blood 
flow, and L inc is an increasing logistic function represented by 
the following equation: 

50 

min -max 

L mc {x, r, max, min) = maxH - : 

1 — min 

1 + -e™ 

max - 1 

55 

wherein max is a maximum value for x, min is a minimum 
value for x, r>0, -oo<x<oo, and min<L. wc . 

12. A medical device for monitoring a circulatory, system, 
the medical device configured to perform a method according 

60 to claim 1 . 

13. A method of monitoring a circulatory system, the 
method comprising: 

measuring a desired circulatory system value; and 
65 comparing said desired circulatory system value to a cor- 
responding circulatory system value calculated using a 
method according to claim 1. 
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14. A non-transitory computer readable medium contain- 
ing computer executable instructions implementing a method 
of modeling a circulatory system, the instructions compris- 
ing: 

a first set of instructions for providing a circulatory system 5 
model including one or more time-dependent pressure 
functions, each of said one or more time-dependent pres- 
sure functions representing a pressure of a portion of the 
circulatory system; 

a second set of instructions for using a logistic function to 10 
represent a regulatory mechanism parameter, said regu- 
latory mechanism parameter representing a regulatory 
mechanism having an impact on circulatory system 
function; 

a third set of instructions for solving one or more equations 1 5 
for said one or more time-dependent pressure functions 
and said logistic function to determine a circulatory 
system dynamic; 

a fourth set of instructions for dividing the circulatory 
system into a plurality of compartments and a heart 20 
pump, said plurality of compartments representing a 
portion of the circulatory system, said heart pump inter- 
acting with at least one of said plurality of compart- 
ments; 

a fifth set of instructions for assigning a time-dependent 25 
pressure function to each of said plurality of compart- 
ments; and 

a sixth set of instructions for using a smooth muscle con- 
traction parameter, said smooth muscle contraction 
parameter representing an autonomic nervous system 30 
effect on smooth muscle contraction in a wall of one of 
said plurality of compartments, wherein said smooth 
muscle contraction parameter is represented according 
to the following equation: 

... 35 

U=C-(P(<j)+P*)) 

wherein Visa time derivative of a volume of said one of 
said plurality of compartments, P is a time derivative 
of the pressure inside said one of said plurality of 
compartments, (j) is a time derivative of said smooth 40 
muscle contraction parameter, P* is a time derivative 
of the pressure outside of said one of said plurality of 
compartments, and C is an active compliance defined 
according to the following equation: 

45 

r yMax e -r{P-(<p+P*)) 

C ~ (1 + e -r{P-{4>+P*))f 

50 

wherein r> 0 , \ rMax is a maximum volume of said one of 
said plurality of compartments, P is the pressure 
inside said one of said plurality of compartments, P* 
is the pressure outside said one of said plurality of 
compartments, and § is said smooth muscle contrac- 55 
tion parameter. 

15. A computer readable medium according to claim 14, 

wherein said fourth set of instructions comprises a seventh set 
of instructions for dividing the circulatory system into an 
arterial compartment and a venous compartment. 60 

16. A computer readable medium according to claim 15, 
further comprising a eighth set of instructions for using a 
logistic function to represent a cerebral blood flow parameter. 

17. A computer readable medium according to claim 16, 
wherein said eighth set of instructions comprises a ninth set of 65 
instructions for representing said cerebral blood flow param- 
eter according to the following logistic function: 


50 

Qi=L inc (PAJ^AV,' 0 - 15 , 1 . 0001 , 0 ) 0 ! 

wherein Q x is said cerebral blood flow parameter, P AV is a 
pressure difference P A -P v , wherein P A is said time-depen- 
dent pressure function for said arterial compartment andP F is 
said time-dependent pressure function for said venous com- 
partment, ^AV is a mean systemic indicative pressure, is a 
mean cerebral blood flow, and L inc is an increasing logistic 
function represented by the following equation: 


min -max 

L mc (x, r, max, min) = maxH 7 — 


wherein max is a maximum value for x, min is a minimum 
value for x, r> 0 , -oo<x<oo, and mm<L inc . 

18. A computer readable medium according to claim 15, 
further comprising a tenth set of instructions for using a 
logistic function to represent an output versus pressure 
parameter, said output versus pressure parameter represent- 
ing an effect of venous pressure on cardiac uptake. 

19. A computer readable medium according to claim 18, 
wherein said tenth set of instructions comprises a eleventh set 
of instructions for representing said output versus pressure 
parameter according to the following logistic function: 

OVP(P v )=L inc (Pjr-P v , 0.5 ,2.5,0) 

wherein OVP is said output versus pressure parameter, P K is 
said time-dependent pressure function for said venous com- 
partment, P v is the mean systemic venous pressure, and L inc is 
an increasing logistic function represented by the following 
equation: 


min — max 

Li nc (x, r, max, min) = maxH : 

1 - min 

1 H -e™ 

max - 1 

wherein max is a maximum value for x, min is a minimum 
value for x, r> 0 , -oo<x<oo, and mm<L irjc . 

20. A computer readable medium according to claim 16, 
wherein said second set of instructions comprises a twelfth 
set of instructions for representing said regulatory mecha- 
nism parameter as an autonomic nervous system effect on 
cardiac output according to the following equation: 

ANSo=L dec (P A -T A , 0.1, 2,0) 

wherein ANS 0 is said regulatory mechanism parameter rep- 
resenting said autonomic nervous system effect on cardiac 
output, P^ is said time-dependent pressure function for said 
arterial compartment, P A is an initial mean arterial pressure, 
and L dec is a decreasing logistic function represented by the 
following equation: 


min - max 

Ldec (A U max, min) = maxH : : 

I — min 

1 + e ~ « 

max - 1 

wherein max is a maximum value for x, min is a minimum 
value for x, r>0, -oo<x<oo, and L^ ec <max. 

21. A computer readable medium according to claim 16, 
wherein said second set of instructions comprises a thirteenth 
set of instructions for representing said regulatory mecha- 
nism parameter as an autonomic nervous system effect on 
non-cerebral blood flow according to the following equation: 
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ANSz=L tn JP A -P A fi3,\.\fil) 
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wherein ANSz is said regulatory mechanism parameter rep- 
resenting said autonomic nervous system effect on non-cere- 
bral blood flow, is said time-dependent pressure function 
for said arterial compartment, is an initial mean arterial 
pressure, and L inc is an increasing logistic function repre- 
sented by the following equation: 


Li n c (x, r, max, min) = max + 


min - max 


1 - min 

1 H -e™ 

max - 1 


wherein max is a maximum value for x, min is a minimum 
value for x, r>0, -oo<x<oo, and mm<L inc . 

22. A computer readable medium according to claim 16, 
wherein said second set of instructions comprises a four- 
teenth set of instructions for representing said regulatory 
mechanism parameter as a central nervous system effects on 
cardiac output according to the following equation: 

CNSo=L dec (Q i -'Q 1 , 0.01,5, 0.9) 

wherein is cerebral blood flow, Q x is mean cerebral blood 
flow, and L dec is a decreasing logistic function represented by 
the following equation: 


min- max 

L dec {x, r, max, min) = maxH : 

1 — min 

1 h 

max - 1 

wherein max is a maximum value for x, min is a minimum 
value for x, r>0, -oo<x<oo, and L dec <max. 

23. A computer readable medium according to claim 16, 
wherein said second set of instructions comprises a fifteenth 
set of instructions for representing said regulatory mecha- 
nism parameter as a central nervous system effects on non- 
cerebral blood flow according to the following equation: 


CNSz^tedQiQi, o.oi, 1.01, o.i) 

wherein Q x is cerebral blood flow, Q x is mean cerebral blood 
flow, and L inc is an increasing logistic function represented by 
the following equation: 


L mc (x, r, max, min) = max + 


min - max 


1 - min 

1 H -e™ 

max - 1 


wherein max is a maximum value for x, min is a minimum 
value for x, r>0, -oo<x<oo, and min<L z „ c . 

24. A method of modeling a pressure and volume relation- 
ship in a compliant vessel using a programmed computer, the 
method comprising: 

defining a first parameter as a change in pressure within the 
vessel using the computer, said change in pressure being 
due to a contraction of smooth muscles of a wall of the 
vessel; 

defining a second parameter as an active compliance for the 
vessel using the computer, said active compliance vary- 
ing with internal pressure, external pressure, and said 
first parameter, wherein said second parameter is 
defined according to the following equation: 


r yMax e -r(P-(<P+P *)) 

Q — 

(1 + e -* P - W +/*))) 2 
5 


wherein C is said second parameter, r>0, y Max is a maxi- 
mum volume of said vessel, P is the pressure inside 
said vessel, P* is the pressure outside said vessel, and 
10 <|) is said first parameter; and 

relating, using the computer, said first parameter and said 
second parameter via a differential equation, wherein 
said relating step comprises relating said first parameter 
and said second parameter according to the following 
15 equation: 

v=o(P-(fc=P*)) 

wherein Visa time derivative of a volume of said vessel, 
P is a time derivative of the pressure inside said vessel, 
20 (j) is a time derivative of said first parameter, P* is a 

time derivative of the pressure outside of said vessel, 
and C is said second parameter. 

25. A method according to claim 24, further comprising 
using said first parameter and said second parameter in mod- 

25 eling a circulatory system. 

26. A non-transitory computer readable storage medium 
containing computer executable instructions implementing a 
method of modeling a pressure and volume relationship in a 
compliant vessel, the instructions comprising: 

30 a first set of instructions for defining a first parameter as a 

change in pressure within the vessel, said change in 
pressure being due to a contraction of smooth muscles of 
a wall of the vessel; 

35 a second set of instructions for defining a second parameter 

as an active compliance for the vessel, said active com- 
pliance varying with internal pressure, external pres- 
sure, and said first parameter, wherein said second set of 
instructions comprises a set of instructions for defining 
4 said second parameter according to the following equa- 

tion: 


r yMax e -riP-(<p+F **)) 

45 C= (1+ c -*p-(^/*))) 2 

wherein C is said second parameter, r>0, y Max is a maxi- 
mum volume of said vessel, P is the pressure inside 
said vessel, P* is the pressure outside said vessel, and 
4> is said first parameter; and 

a third set of instructions for relating said first parameter 
and said second parameter via a differential equation, 
wherein said third set of instructions comprises a set of 
55 instructions for relating said first parameter and said 
second parameter according to the following equation: 

r=ciP-(faP*)) 

wherein Visa time derivative of a volume of said vessel, 
60 Pisa time derivative of the pressure inside said vessel , 

(j) is a time derivative of said first parameter, P* is a 
time derivative of the pressure outside of said vessel, 
and C is said second parameter. 

27. A computer readable medium according to claim 26, 
65 further comprising a sixth set of instructions for using said 
first parameter and said second parameter in modeling a cir- 
culatory system. 
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28. A method of modeling a circulatory system using a 
programmed computer, the method comprising: 

providing a circulatory system model in the computer, the 
model including one or more time-dependent pressure 
functions, each of said one or more time-dependent pres- 5 
sure functions representing a pressure in a portion of the 
circulatory system; 

using a logistic function in the computer to represent a 
regulatory mechanism parameter, said regulatory 
mechani sm parameter representing a regulatory mecha- 1 0 
nism having an impact on circulatory system function, 
said using a logistic function including using a logistic 
function to represent a cerebral blood flow parameter 
according to the following logistic function: 

21=^(^^0-15,1.0001, 0)-5i 

wherein Q x is said cerebral blood flow parameter, P^is 
a pressure difference P A -P V , wherein is said time- 
dependent pressure function for said arterial compart- 
ment and P K is said time-dependent pressure function 20 
for said venous compartment, P^is a mean systemic 
indicative pressure, Q x is a mean cerebral blood flow, 
and L inc is an increasing logistic function represented 
by the following equation: 


min - max 

Lmc (u r, max, min) = max H ; : 

1 - min 

1 H -e™ 

max - 1 

30 

wherein max is a maximum value for x, min is a mini- 
mum value for x, r>0, -oo<x<oo, and min<L zwc ; 

dividing the circulatory system into a plurality of compart- 
ments and a heart pump, said plurality of compartments 
representing a portion of the circulatory system, said 35 
heart pump interacting with at least one of said plurality 
of compartments, wherein said dividing step divides 
said circulatory system into a set of compartments com- 
prising of an arterial compartment and a venous com- 
partment, and a heart pump; 

assigning, using the computer, one of said one or more 
time-dependent pressure functions to each of said plu- 
rality of compartments; 

using a smooth muscle contraction parameter, said smooth 
muscle contraction parameter representing an auto- 
nomic nervous system effect on smooth muscle contrac- 
tion in a wall of one of said plurality of compartments; 
and 

solving, using the computer, one or more equations for said 
one or more time-dependent pressure functions and said 
logistic function to determine a circulatory system 
dynamic. 

29. A method of modeling a circulatory system using a 
programmed computer, the method comprising: 

providing a circulatory system model in the computer, the 
model including one or more time-dependent pressure 
functions, each of said one or more time-dependent pres- 
sure functions representing a pressure in a portion of the 
circulatory system; 

using a logistic function in the computer to represent a 
regulatory mechanism parameter, said regulatory 
mechanism parameter representing a regulatory mecha- 
nism having an impact on circulatory system function, 
said using a logistic function including using a logistic 
function to represent an output versus pressure param- 
eter, said output versus pressure parameter representing 


54 

an effect of venous pressure on cardiac uptake according 
to the following logistic function: 

OVP{PJ)=L inc {Pv-P v , 0.5, 2.5,0) 

wherein OVP is said output versus pressure parameter, 
P K is said time-dependent pressure function for said 
venous compartment, P K is the mean systemic venous 
pressure, and L inc is an increasing logistic function 
represented by the following equation: 


min - max 

L'mc (x, r, max, min) = maxH : 

1 — min 

1 + -e™ 

max - 1 

wherein max is a maximum value for x, min is a mini- 
mum value for x, r>0, -oo<x<oo ? and min<L zwc ; 

dividing the circulatory system into a plurality of compart- 
ments and a heart pump, said plurality of compartments 
representing a portion of the circulatory system, said 
heart pump interacting with at least one of said plurality 
of compartments, wherein said dividing step divides 
said circulatory system into a set of compartments com- 
prising of an arterial compartment and a venous com- 
partment, and a heart pump; 

assigning, using the computer, one of said one or more 
time-dependent pressure functions to each of said plu- 
rality of compartments; 

using a smooth muscle contraction parameter, said smooth 
muscle contraction parameter representing an auto- 
nomic nervous system effect on smooth muscle contrac- 
tion in a wall of one of said plurality of compartments; 
and 

solving, using the computer, one or more equations for said 
one or more time-dependent pressure functions and said 
logistic function to determine a circulatory system 
dynamic. 

30. A method of modeling a circulatory system using a 
programmed computer, the method comprising: 

providing a circulatory system model in the computer, the 
model including one or more time-dependent pressure 
functions, each of said one or more time-dependent pres- 
sure functions representing a pressure in a portion of the 
circulatory system; 

using a logistic function in the computer to represent a 
regulatory mechanism parameter, said regulatory 
mechanism parameter representing a regulatory mecha- 
nism having an impact on circulatory system function, 
wherein said regulatory mechanism parameter repre- 
sents an autonomic nervous system effect on cardiac 
output according to the following equation: 

ANSo=L dec (P A -P A , 0 . 1 , 2 , 0 ) 

wherein ANS 0 is said regulatory mechanism parameter 
representing said autonomic nervous system effect on 
cardiac output, P A is said time-dependent pressure 
function for said arterial compartment, P A is an initial 
mean arterial pressure, and L dec is a decreasing logis- 
tic function represented by the following equation: 


min - max 

L dec (x, r, max, min) = maxH : 

1 - mm 

1 H -e ™ 

max - 1 



US 7,794,230 B2 


55 

wherein max is a maximum value for x, min is a mini- 
mum value for x, r>0, -oo<x<oo ; and L^ ec <max; 
dividing the circulatory system into a plurality of compart- 
ments and a heart pump, said plurality of compartments 
representing a portion of the circulatory system, said 5 
heart pump interacting with at least one of said plurality 
of compartments, wherein said dividing step divides 
said circulatory system into a set of compartments com- 
prising of an arterial compartment and a venous com- 
partment, and a heart pump; to 

assigning, using the computer, one of said one or more 
time-dependent pressure functions to each of said plu- 
rality of compartments; 

using a smooth muscle contraction parameter, said smooth 
muscle contraction parameter representing an auto- 15 
nomic nervous system effect on smooth muscle contrac- 
tion in a wall of one of said plurality of compartments; 
and 

solving, using the computer, one or more equations for said 
one or more time-dependent pressure functions and said 
logistic function to determine a circulatory system 
dynamic. 

31. A method of modeling a circulatory system using a 
programmed computer, the method comprising: 25 

providing a circulatory system model in the computer, the 
model including one or more time -dependent pressure 
functions, each of said one or more time-dependent pres- 
sure functions representing a pressure in a portion of the 
circulatory system; 30 

using a logistic function in the computer to represent a 
regulatory mechanism parameter, said regulatory 
mechanism parameter representing a regulatory mecha- 
nism having an impact on circulatory system function, 
wherein said regulatory mechanism parameter repre- 35 
sents an autonomic nervous system effect on non-cere- 
bral blood flow according to the following equation: 

ANSz=L inc (P A -P A , 0.3,1 . 1 ,0.7) 

wherein ANSz is said regulatory mechanism parameter 40 
representing said autonomic nervous system effect on 
non-cerebral blood flow, p* is said time-dependent pres- 
sure function for said arterial compartment, P A is an 
initial mean arterial pressure, and L inc is an increasing 
logistic function represented by the following equation: 45 


min - max 

Line U, U max, min) = max H : — 

1 - min 


wherein max is a maximum value for x, min is a mini- 
mum value for x, r>0, -oo<x<oo, and min<L z „ c ; 
dividing the circulatory system into a plurality of compart- 55 
ments and a heart pump, said plurality of compartments 
representing a portion of the circulatory system, said 
heart pump interacting with at least one of said plurality 
of compartments, wherein said dividing step divides 
said circulatory system into a set of compartments com- 60 
prising of an arterial compartment and a venous com- 
partment, and a heart pump; 

assigning, using the computer, one of said one or more 
time-dependent pressure functions to each of said plu- 
rality of compartments; 65 

using a smooth muscle contraction parameter, said smooth 
muscle contraction parameter representing an auto- 
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nomic nervous system effect on smooth muscle contrac- 
tion in a wall of one of said plurality of compartments; 
and 

solving, using the computer, one or more equations for said 
one or more time-dependent pressure functions and said 
logistic function to determine a circulatory system 
dynamic. 

32. A method of modeling a circulatory system using a 
programmed computer, the method comprising: 

providing a circulatory system model in the computer, the 
model including one or more time-dependent pressure 
functions, each of said one or more time-dependent pres- 
sure functions representing a pressure in a portion of the 
circulatory system; 

using a logistic function in the computer to represent a 
regulatory mechanism parameter, said regulatory 
mechanism parameter representing a regulatory mecha- 
nism having an impact on circulatory system function, 
wherein said regulatory mechanism parameter repre- 
sents a central nervous system effect on cardiac output 
according to the following equation: 

CNSo=L dec (Qi~Qi, 0.01, 5, 0.9) 

wherein Q x is cerebral blood flow, Q x is mean cerebral 
blood flow, and L dec is a decreasing logistic function 
represented by the following equation: 


min - max 

Ldecix, r, max, min) = maxH ; : 

1 - min 

1 H -e rx 

max - 1 

wherein max is a maximum value for x, min is a mini- 
mum value for x, r>0, -oo<x<oo ? and L rfec <max; 

dividing the circulatory system into a plurality of compart- 
ments and a heart pump, said plurality of compartments 
representing a portion of the circulatory system, said 
heart pump interacting with at least one of said plurality 
of compartments, wherein said dividing step divides 
said circulatory system into a set of compartments com- 
prising of an arterial compartment and a venous com- 
partment, and a heart pump; 

assigning, using the computer, one of said one or more 
time-dependent pressure functions to each of said plu- 
rality of compartments; 

using a smooth muscle contraction parameter, said smooth 
muscle contraction parameter representing an auto- 
nomic nervous system effect on smooth muscle contrac- 
tion in a wall of one of said plurality of compartments; 
and 

solving, using the computer, one or more equations for said 
one or more time-dependent pressure functions and said 
logistic function to determine a circulatory system 
dynamic. 

33. A method of modeling a circulatory system using a 
programmed computer, the method comprising: 

providing a circulatory system model in the computer, the 
model including one or more time-dependent pressure 
functions, each of said one or more time-dependent pres- 
sure functions representing a pressure in a portion of the 
circulatory system; 

using a logistic function in the computer to represent a 
regulatory mechanism parameter, said regulatory 
mechanism parameter representing a regulatory mecha- 
nism having an impact on circulatory system function, 
wherein said regulatory mechanism parameter repre- 
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sents a central nervous system effects on non-cerebral 
blood flow according to the following equation: 

CNSz=L^(0i-2 l5 0.01, 1.01 ,0.1) 

wherein Q x is cerebral blood flow, Q x is mean cerebral 
blood flow, and h inc is an increasing logistic function 
represented by the following equation: 


Line (jf) r, max, min) = max + 


min- max 


1 4- 


1 - min 

r e ' 

max - 1 


10 


wherein max is a maximum value for x, min is a mini- 15 
mum value for x, r>0, -oo<x<oo, and min<L IWC ; 

dividing the circulatory system into a plurality of compart- 
ments and a heart pump, said plurality of compartments 
representing a portion of the circulatory system, said 20 
heart pump interacting with at least one of said plurality 
of compartments, wherein said dividing step divides 
said circulatory system into a set of compartments com- 
prising of an arterial compartment and a venous com- 
partment, and a heart pump; 

assigning, using the computer, one of said one or more 
time-dependent pressure functions to each of said plu- 
rality of compartments; 

using a smooth muscle contraction parameter, said smooth 30 
muscle contraction parameter representing an auto- 
nomic nervous system effect on smooth muscle contrac- 
tion in a wall of one of said plurality of compartments; 
and 

solving, using the computer, one or more equations for said 
one or more time-dependent pressure functions and said 
logistic function to determine a circulatory system 
dynamic. 

34 . A non-transitory computer readable medium contain- 40 
ing computer executable instructions implementing a method 
of modeling a circulatory system, the instructions compris- 
ing: 

a set of instructions for providing a circulatory system 
model including one or more time-dependent pressure 45 
functions, each of said one or more time-dependent pres- 
sure functions representing a pressure of a portion of the 
circulatory system; 

a set of instructions for using a logistic function to repre- 50 
sent a regulatory mechanism parameter, said regulatory 
mechanism parameter representing a regulatory mecha- 
nism having an impact on circulatory system function; 

a set of instructions for using a logistic function to repre- 
sent a cerebral blood flow parameter according to the 55 
following logistic function: 

0.15,1.0001,0)-^ 

wherein Q x is said cerebral blood flow parameter, P^is 6Q 
a pressure difference V A -P v , wherein P A is said time- 
dependent pressure function for said arterial compart- 
ment and P F is said time-dependent pressure function 
for said venous compartment, P AV is a mean systemic 
indicative pressure, Q : is a mean cerebral blood flow, 65 
and L inc is an increasing logistic function represented 
by the following equation: 


L mc {x, r, max, min) = max + 


min-max 


1 4- 


1 - min 

e rx 

max - 1 


wherein max is a maximum value for x, min is a mini- 
mum value for x, r>0, -oo<x<oo, and min<L zwc a set of 
instructions for solving one or more equations for said 
one or more time-dependent pressure functions and 
said logistic function to determine a circulatory sys- 
tem dynamic; 

a set of instructions for dividing the circulatory system into 
a plurality of compartments and a heart pump, said plu- 
rality of compartments representing a portion of the 
circulatory system, said heart pump interacting with at 
least one of said plurality of compartments, wherein the 
plurality of compartments includes an arterial compart- 
ment and a venous compartment; 

a set of instructions for assigning a time-dependent pres- 
sure function to each of said plurality of compartments; 
and 

a set of instructions for using a smooth muscle contraction 
parameter, said smooth muscle contraction parameter 
representing an autonomic nervous system effect on 
smooth muscle contraction in a wall of one of said plu- 
rality of compartments. 

35. A non-transitory computer readable medium contain- 
ing computer executable instructions implementing a method 
of modeling a circulatory system, the instructions compris- 
ing: 

a set of instructions for providing a circulatory system 
model including one or more time-dependent pressure 
functions, each of said one or more time-dependent pres- 
sure functions representing a pressure of a portion of the 
circulatory system; 

a set of instructions for using a logistic function to repre- 
sent a regulatory mechanism parameter, said regulatory 
mechanism parameter representing a regulatory mecha- 
nism having an impact on circulatory system function; 

a set of instructions for using a logistic function to repre- 
sent an output versus pressure parameter that represents 
an effect of venous pressure on cardiac uptake wherein 
the representing the output versus pressure parameter 
includes the following logistic function: 

OYY{P v )=L inc {P r -T v ,0.5, 2.5,0) 

wherein OVP is said output versus pressure parameter, 
P v is said time-dependent pressure function for said 
venous compartment, P F is the mean systemic venous 
pressure, and L zwc is an increasing logistic function 
represented by the following equation: 


L inc (x, r, max, min) = max 4- 


min-max 


1 - min 

1 H -e™ 

max - 1 


wherein max is a maximum value for x, min is a mini- 
mum value for x, r>0, -oo<x<oo 5 and min<L zwc a set of 
instructions for solving one or more equations for said 
one or more time-dependent pressure functions and 
said logistic function to determine a circulatory sys- 
tem dynamic; 

a set of instructions for dividing the circulatory system into 
a plurality of compartments and a heart pump, said plu- 
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rality of compartments representing a portion of the 
circulatory system, said heart pump interacting with at 
least one of said plurality of compartments, wherein the 
plurality of compartments includes an arterial compart- 
ment and a venous compartment; 

a set of instructions for assigning a time-dependent pres- 
sure function to each of said plurality of compartments; 
and 

a set of instructions for using a smooth muscle contraction 
parameter, said smooth muscle contraction parameter 
representing an autonomic nervous system effect on 
smooth muscle contraction in a wall of one of said plu- 
rality of compartments. 

36. A non-transitory computer readable medium contain- 
ing computer executable instructions implementing a method 
of modeling a circulatory system, the instructions compris- 
ing: 

a set of instructions for providing a circulatory system 
model including one or more time-dependent pressure 
functions, each of said one or more time-dependent pres- 
sure functions representing a pressure of a portion of the 
circulatory system; 

a set of instructions for using a logistic function to repre- 
sent a regulatory mechanism parameter, said regulatory 
mechanism parameter representing a regulatory mecha- 
nism having an impact on circulatory system function, 
said regulatory mechanism parameter representing an 
autonomic nervous system effect on cardiac output 
according to the following equation: 

ANSo=L dec (P A -P A ,0A,2,0) 

wherein ANS 0 is said regulatory mechanism parameter 
representing said autonomic nervous system effect on 
cardiac output, P A is said time-dependent pressure 
function for said arterial compartment, P A is an initial 
mean arterial pressure, and h dec is a decreasing logis- 
tic function represented by the following equation: 

min -max 

Ld ec (x, r, max, min) = max H : 

1 - min 

1 H -e rx 

max - 1 

wherein max is a maximum value for x, min is a mini- 
mum value for x, r>0, -oo<x<oo 5 and L rfec <max; 

a set of instructions for solving one or more equations for 
said one or more time-dependent pressure functions and 
said logistic function to determine a circulatory system 
dynamic; 

a set of instructions for dividing the circulatory system into 
a plurality of compartments and a heart pump, said plu- 
rality of compartments representing a portion of the 
circulatory system, said heart pump interacting with at 
least one of said plurality of compartments, wherein the 
plurality of compartments includes an arterial compart- 
ment and a venous compartment; 

a set of instructions for assigning a time-dependent pres- 
sure function to each of said plurality of compartments; 
and 

a set of instructions for using a smooth muscle contraction 
parameter, said smooth muscle contraction parameter 
representing an autonomic nervous system effect on 
smooth muscle contraction in a wall of one of said plu- 
rality of compartments. 
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37. A non-transitory computer readable medium contain- 
ing computer executable instructions implementing a method 
of modeling a circulatory system, the instructions compris- 
ing: 

5 a set of instructions for providing a circulatory system 
model including one or more time-dependent pressure 
functions, each of said one or more time-dependent pres- 
sure functions representing a pressure of a portion of the 
circulatory system; 

to a set of instructions for using a logistic function to repre- 
sent a regulatory mechanism parameter, said regulatory 
mechanism parameter representing a regulatory mecha- 
nism having an impact on circulatory system function, 
said regulatory mechanism parameter representing an 
15 autonomic nervous system effect on non-cerebral blood 
flow according to the following equation: 

ANSz=L inc {P A -PA$3 ,1 ■ 1 , 0 . 7 ) 

wherein ANSz is said regulatory mechanism parameter 
20 representing said autonomic nervous system effect on 

non-cerebral blood flow, P A is said time-dependent 
pressure function for said arterial compartment, P A is 
an initial mean arterial pressure, and L inc is an increas- 
ing logistic function represented by the following 
25 equation: 

min - max 

L mc (x, r, max, min) = maxH ; : 

1 - min 
1 H -e™ 

30 max - 1 

wherein max is a maximum value for x, min is a mini- 
mum value for x, r>0, -oo<x<oo 5 and min<L /wc . 

a set of instructions for solving one or more equations for 
said one or more time-Dependent pressure functions and 
said logistic function to determine a circulatory system 
dynamic; 

a set of instructions for dividing the circulatory system into 
a plurality of compartments and a heart pump, said plu- 
rality of compartments representing a portion of the 
circulatory system, said heart pump interacting with at 
least one of said plurality of compartments, wherein the 
plurality of compartments includes an arterial compart - 
ment and a venous compartment; 

a set of instructions for assigning a time-dependent pres- 
sure function to each of said plurality of compartments; 
and 

a set of instructions for using a smooth muscle contraction 
50 parameter, said smooth muscle contraction parameter 
representing an autonomic nervous system effect on 
smooth muscle contraction in a wall of one of said plu- 
rality of compartments. 

38. A non-transitory computer readable medium contain- 
55 ing computer executable instructions implementing a method 

of modeling a circulatory system, the instructions compris- 
ing: 

a set of instructions for providing a circulatory system 
model including one or more time-dependent pressure 
60 functions, each of said one or more time-dependent pres- 

sure functions representing a pressure of a portion of the 
circulatory system; 

a set of instructions for using a logistic function to repre- 
sent a regulatory mechanism parameter, said regulatory 
65 mechanism parameter representing a regulatory mecha- 
nism having an impact on circulatory system function, 
said regulatory mechanism parameter representing an a 
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central nervous system effects on cardiac output accord- a set of instructions for using a logistic function to repre- 

ing to the following equation: sent a regulatory mechanism parameter, said regulatory 

CNSo^L dec {Q ) -'Q^ o.oi ,5,0.9) mechanism parameter representing a regulatory mecha- 


wherein Q x is cerebral blood flow, Q x is mean cerebral 
blood flow, and L dec is a decreasing logistic function 
represented by the following equation: 


min -max 

Ld ec (x, r, max, min) = max H ; 

1 — min 


wherein max is a maximum value for x, min is a mini- 
mum value for x, r>0, -oo<x<oo, and L dec <max; 

a set of instructions for solving one or more equations for 
said one or more time-dependent pressure functions and 
said logistic function to determine a circulatory system 
dynamic; 

a set of instructions for dividing the circulatory system into 
a plurality of compartments and a heart pump, said plu- 
rality of compartments representing a portion of the 
circulatory system, said heart pump interacting with at 
least one of said plurality of compartments, wherein the 
plurality of compartments includes an arterial compart- 
ment and a venous compartment; 

a set of instructions for assigning a time-dependent pres- 
sure function to each of said plurality of compartments; 
and 

a set of instructions for using a smooth muscle contraction 
parameter, said smooth muscle contraction parameter 
representing an autonomic nervous system effect on 
smooth muscle contraction in a wall of one of said plu- 
rality of compartments. 

39. A non-transitory computer readable medium contain- 
ing computer executable instructions implementing a method 
of modeling a circulatory system, the instructions compris- 
ing: 

a set of instructions for providing a circulatory system 
model including one or more time-dependent pressure 
functions, each of said one or more time-dependent pres- 
sure functions representing a pressure of a portion of the 
circulatory system; 


nism having an impact on circulatory system function, 
said regulatory mechanism parameter representing an a 
central nervous system effects on non-cerebral blood 
flow according to the following equation: 

CNSz=L inc (Qi~Qi, 0.01,1.01,0.1) 

wherein Q x is cerebral blood flow, Q x is mean cerebral 
blood flow, and L inc is an increasing logistic function 
represented by the following equation: 


nun - max 

L inc (x, r, max, min) = maxH : 

1 — min 


20 wherein max is a maximum value for x, min is a mini- 
mum value for x, r>0, -oo<x<oo, and min<L /wc ; 
a set of instructions for solving one or more equations for 
said one or more time-dependent pressure functions and 
said logistic function to determine a circulatory system 
25 dynamic; 

a set of instructions for dividing the circulatory system into 
a plurality of compartments and a heart pump, said plu- 
rality of compartments representing a portion of the 
circulatory system, said heart pump interacting with at 
least one of said plurality of compartments, wherein the 
plurality of compartments includes an arterial compart- 
ment and a venous compartment; 
a set of instructions for assigning a time-dependent pres- 
„ sure function to each of said plurality of compartments; 

and 

a set of instructions for using a smooth muscle contraction 
parameter, said smooth muscle contraction parameter 
representing an autonomic nervous system effect on 
40 smooth muscle contraction in a wall of one of said plu- 
rality of compartment. 



